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Introduction 



An extensive search for top quarks has been performed at electron-positron and hadron 
coUiders for more than a decade. First evidence for top quark production in proton- 
antiproton colhsions has been announced by the CDF collaboration in the spring of 1994. 
After collecting more luminosity subsequently both the CDF and the DO experiments 
presented the definite analysis which demonstrated not only the existence of the anti- 
cipated quark but at the same time also provided a kinematic determination of the top 
quark mass around 180 GeV and a production cross section consistent with the QCD pre- 
dictions. The mass value is in perfect agreement with the indirect mass determinations 
based on precision measurements [2-7] of the electroweak parameters in e^e~ annihilation 
and in lepton-nucleon scattering. Exploiting the quadratic top mass dependence of radia- 
tive corrections an indirect mass measurement of 180 GeV with a present uncertainty of 
roughly 20 GeV has been achieved. 

The top quark completes the fermionic spectrum of the Standard Model. All its properties 
are uniquely fixed after the mass has been determined. However, as a consequence of its 
large mass and decay rate it will behave markedly different compared to the remaining 
five lighter quarks. 

It is not just the obvious aim for completion which raises the interest in the top quark. 
With its mass comparable to the scale of the electroweak symmetry breaking it is plausible 
that top quark studies could provide an important clue for the understanding of fermion 
mass generation and the pattern of Yukawa couplings. In fact, it has been suggested 
that a top quark condensate could even be responsible for the mechanism of spontaneous 
symmetry breaking p[. 

These lectures will be mainly concerned with top quark phenomenology within the frame- 
work of the Standard Model (SM). The precise understanding of its production and decay 
constitutes the basis of any search for deviations or physics beyond the SM. 

The properties of the top quark will be covered in chapter 1. Direct and indirect determi- 
nations of its decay rates, decay distributions including QCD and electroweak corrections 
and decay modes predicted in supersymmetric extensions will be discussed. Top quark 
production at hadron colliders will be the subject of chapter 2. The production cross sec- 
tion and momentum distributions are important ingredients in any of the present analysis. 
An alternative reaction, namely top quark production through H^-6-fusion allows to de- 
termine the W-b-i coupling and thus indirectly the top quark decay rate. 

Perspectives for top studies at a future e^e~collider will be presented in chapter 3. An 
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accurate determination of the top quark mass and its width to better than 1 GeV with 
a relative accuracy of about 10% seems feasible, and the electroweak couplings of the 
top quark could be precisely measured with the help of polarized beams. Of particular 
interest is the interplay between the large top quark decay rate and the binding through 
the QCD potential which will be also covered in chapter 3. 
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Chapter 1 

The Profile of the Top Quark 



Hadron collider experiments at the TEVATRON have firmly established the existence 
of the top quark and already provide a fairly accurate determination of its mass. The 
couplings of the top quark to the gauge bosons are uniquely fixed by the SM. Thus all its 
properties — its production cross section and its decay rate and distributions — can be 
predicted unambiguously. 

The study of real top quarks at high energy colliders, in particular the observation of a 
peak in the invariant mass of its decay products, is certainly the most impressive proof 
of existence. Nevertheless, the indirect evidence for a top quark and the determination of 
its mass is not only of historical interest. The arguments which anticipated the existence 
of the top quark and its mass around 180 GeV illustrate the rigid structure of the SM, 



its selfconsistency and beauty. They will be presented in section pTl . 

These theoretical arguments have inspired the experimental searches. The upper limit 
on the top mass around 200 GeV deduced already relatively early from electroweak pre- 
cision studies has provided encouragement that energies of present coUiders were suited 
to complete this enterprise. The agreement between the most recent indirect mass deter- 
minations both through radiative corrections and through direct observation strengthens 
the present belief into the quantum field aspect of the theory. It furthermore justifies the 
corresponding line of reasoning concerning the search for the last remaining ingredient of 



the SM, the Higgs boson. Section |L1| of this first chapter will, with this motivation in 
mind, be devoted to a discussion of the indirect information on the top quark, its existence 
and its mass. Top decays, including aspects of radiative corrections, polarisation effects 



and decays induced by physics beyond the SM will be covered in section Q]2 



1.1 Indirect information 
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1.1.1 Indirect evidence for the top quark 

Several experimental results already prior to its discovery did provide strong evidence 
that the fermion spectrum of the Standard Model 
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does include the top quark, imprinting the same multiplet structure on the third family 
as the first two families. The evidence is based on theoretical selfconsistency (absence of 
anomalies), the absence of flavour changing neutral currents (FCNC) and measurements 
of the weak isospin of the b quark which has been proved to be non-zero, — —1/2, thus 
demanding an /a = +1/2 partner in this isospin multiplet. 



Absence of triangle anomalies 

A compelling argument for the existence of top quarks follows from a theoretical consis- 
tency requirement. The renormalizability of the Standard Model demands the absence of 
triangle anomalies. Triangular fermion loops built-up by an axialvcctor charge I^a = ~hL 
combined with two electric vector charges Q would spoil the renormalizability of the gauge 
theory. Since the anomalies do not depend on the masses of the fermions circulating in 
the loops, it is sufficient to demand that the sum 




of all contributions be zero. Such a requirement can be translated into a condition on the 
electric charges of all the left-handed fermions 



E^3 = 0- 

L 



(1.1) 



This condition is met in a complete standard family in which the electric charges of the 
lepton plus those of all color components of the up and down quarks add up to zero. 



^g = -i + 3x 

L 



0. 



If the top quark were absent from the third family, the condition would be violated and 
the Standard Model would be theoretically inconsistent. 
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Absence of FCNC decays 

Mixing between quarks which belong to different isospin multiplets 



b'r 



b'r 



sl cos-d' + hi sin-i?' 
—sism'd' + 1)1 cos "d' 



generates non-diagonal neutral current couplings, i.e. the breaking of the GIM mechanism 



< /3 > 



'2 lCL,CLj 



2 l^L'^L 



^ (cl, cl) ^ (sl, Sl) cos^ ~\ bi) sin^ ^' 

~ sin^' cos^' ({sl, bi) + sl)) • 



The non-diagonal current induces flavor-changing neutral lepton pair decays b—>- s + 
which have been estimated to be a substantial fraction of all semileptonic B meson decays. 
The relative strenth of neutral versus charged current induced rate is essentially given by 



NC 
^CC 



1 fM'^Y{v! + al){v',+al) 



2 VM| 



0.06. 



;i.2) 



Taking the proper momentum dependence of the matrix element and the phase space into 
account one finds 

BR{B^l+i^iX) - ^ ' 

This ratio is four orders of magnitude larger than a bound set by the UAl Collaboration 

BO 

BR {B^ fi+fi-X) 5.0 X 10-5 

BR {B-^ fiu^X) ^ 0.103 ±0.005' ^ ' 

so that the working hypothesis of an isosinglet b quark is clearly ruled out experimentally 
also by this method. 



Partial width r(Z— > bb) and forward-backward asymmetry of b quarks 

The Z boson couples to quarks through vector and axial-vector charges with the well- 
known strength 




7/. K - ag75] 
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Depending on the isospin assignment of righthanded and lefthanded quark fields these 
charges are defined as 



Vg = 2(/f + If)- 4e,sin2,9i^ (1.5) 
S = 2(/f-/f) (1.6) 

For the present application the Born approximation in the massless limit provides an 
adequate representation of the partial Z decay rate 

rsiZ^bb)^^(3{vl + al). (1.7) 

In the Standard Model 2J^^ = and 21^^ = ±1 for up/down quarks respectively. 

The ratio between the predictions in the context of a topless model and the SM amounts 
to 

ptopless ^ (4g,sin^gH^)^ 

rsM (i + 4g,sin2M' + i -^'/'^ ^'■''> 

whereas theory and LEP experiments are well consistent 

T{Z ^ hh) _ { 0.2155 ± 0.0004 theory [3] , . 

T{Z had) ~ \ 0.2219 ± 0.0017 experiment § ^ ' ' 

ruling out the / = assignement for the 6-quark. The forward-backward asymmetry at 
the Z resonance 

A™ = -^AM (1.10) 

with 

is sensitive toward the relative size of vector and axial h quark couplings. Up to a sign, 
the first of these factors, Ae ~ 0.15, can be interpreted as the degree of longitudinal Z 
polarisation, Pz = —A^, which is induced by the electron coupling even for unpolarized 
beams. For longitudinally polarized beams it can be replaced by unity. The second factor 
represents essentially the analyzing power of b quarks. With a predicted value of 0.93 it 
is close to its maximum in the SM. In fact, this high analysing power is the reason for 
the large sensitivity of A^^ toward sin^ 9w ||12[- ^ fictitious topless model Ab is zero. 



The most recent experimental results from LEP and SLC are displayed in Fig. |1.1| . 

A remaining sign ambiguity is finally resolved by the interference between NC and elec- 
tromagnetic amplitude. It leads to a forward backward asymmetry at low energies 
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World Measurements 
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Figure 1.1: Experimental results for the asymmetry parameter Ai, ||TB 
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which has been studied in particular at PEP, PETRA, and most recently with highest 



precision at Tristan at a cm energy of 58 GeV [Q. Using the data available shortly after 
the turnon of LEP and combining r^j, A^, and ApB 



exp 



exp 



-0.490t°:°l^ 
-0.028±0.056 



I^ib) = -1/2 







has been obtained already some time ago |jT5[ . As shown in Fig. |1.2| all measurements are 
nicely consistent with the predictions of the SM[|. The isospin assignment of the Standard 
Model is thus well confirmed. 



Figure 1.2: The weak isospins I^{b) and I-^{b) of the left- and right-handed b quark 
components, extracted from the data on T{Z^ bb) and Apsib) at LEP, and PETRA/PEP 
and TRISTAN, Ref.0. 



For a discussion of the most recent results for Rb, however, see 
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1.1.2 Mass limits and indirect mass determinations 



Theoretical constraints 



Present theoretical analyses of the Standard Model are based almost exclusively on per- 
turbation theory. If this method is assumed to apply also to the top-quark sector, in 
particular when linked to the Higgs sector, the top mass must be bounded as the strength 
of the Yukawa- Higgs-top coupling is determined by this parameter. The following consis- 
tency conditions must be met: 

Perturbative Yukawa coupling gy (ttH) 



Defining the Yukawa coupling in the Standard Model through 



C 



Y 



9y 



V2 



R 



h.c. 



;i.i3) 



the coupling constant gy is related to the top mass by 

gy {UH) = mt\j2V2GF- 



;i.i4) 



Demanding the effective expansion parameter gy/^Tc to be smaller than 1, the top mass 
is bounded to 



m < 



V27V 



620GeV. 



;i.i5) 



For a top mass of 180 GeV the coupling gy/^TT 0.085 is comfortably small so that 
perturbation theory can safely be applied in this region. 

Unitarity bound 

At asymptotic energies the amplitude of the zeroth partial wave for elastic tt scattering 



in the color singlet same-helicity channel 

ao (tt— >■ tt) = 



57r 



;i.i6) 



TC 



grows quadratically with the top mass. Unitarity however demands this real amplitude 
to be bounded by |Reao| < 1/2. This condition translates to 



mt < 



500 GeV. 



;i.i7) 



The bound improves by taking into account the running of the Yukawa coupling [T^. 
These arguments are equally applicable for any additional species of chiral fermions with 
mass induced via spontaneous symmetry breaking. 
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Stability of the Higgs system: top-Higgs bound 



The quartic couphng A in the effective Higgs potential 

v = ^^'\<P\' + \\4>\' 

depends on the scale at which the system is interacting. The running of A is induced by 
higher-order loops built-up by the Higgs particles themselves, the vector bosons and the 
fermions in the Standard Model |T^. For moderate values of the top mass, nit < 77 



GeV, these radiative corrections would have generated a lower bound on the Higgs mass 
of 7 GeV. With the present experimental lower limit itih > 65 GeV and the top quark 
mass determined around 180 GeV this bound is of no practical relevance any more. At 
high energies the radiative corrections make A rise up to the Landau pole at the cut- 
off parameter A beyond which the Standard Model in the present formulation cannot be 
continued ["triviality bound", as this bound could formally be misinterpreted as requiring 
the low energy coupling to vanish]. If for a fixed Higgs mass the top mass is increased, 
the top loop radiative corrections lead to negative values of the quartic coupling A 

■^T~~ = ^ f-^^ ~ ^9y + gauge couplings] (1.18) 
a logs Sn^ L J 

Since the potential is unbounded from below in this case, the Higgs system becomes 
instable. Thus the stability requirement defines an upper value of the top mass for a 
given Higgs mass rriH and a cut-off scale A. The result of such an analysis is presented in 
Fig. |1.3| . Depending on the cut-off scale A where new physics may set in, the top mass is 
bounded to nit < 200 GeV if A exceeds the Planck scale but it rises up to 400 to 500 GeV 
if the cut-off is reached at a level of 1 TeV and below. The estimates are similar to the 
unitarity analysis in the preceeding subsection. Lattice simulations of the Yukawa model 
have arrived at qualitatively similar results (see e.g. 



These theoretical analyses have shown that for the top mass around 180 GeV the Standard 
Model may be valid up to a cut-off at the Planck scale. [The hierarchy problem, that is 
not touched in the present discussion, may enforce nevertheless new additional physical 
phenomena already in the TeV range.] 

In the context of the SM the top Yukawa coupling is simply present as a free parameter. 
In the minimal supersymmetric model, however, the picture is changed completely. A 
large Yukawa coupling may play the role of a driving term for the spontaneous breaking 



of SU{2) X U{1), as discussed in |2y] and in fact the large mass of the top quark had been 



predicted on the basis of these arguments prior to its experimental discovery. 



Mass estimates from radiative corrections 

First indications of a high top quark mass were derived from the rapid {B — B) oscillations 



observed by ARGUS |21|]. However, due to the uncertainties of the CKM matrix element 
Vtd and of the (M) wave function, not more than qualitative conclusions can be drawn 
from such an analysis as the oscillation frequency Am ~ \Vtdf fs^nt depends on three 
[unknown] parameters. 
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Figure 1.3: Bounds on the Higgs and top masses following from triviality of the Higgs's 
quartic self-coupling and the stability of the Higgs system; from [|l^ . 



The analysis of the radiative corrections to high precision electroweak observables is much 
more advanced [2-7]. Since Higgs mass effects are weak as a result of the screening theorem 



2^ , the top mass is the dominant unknown in the framework of the Standard Model. 
Combining the high precision measurements of the Z mass with sin^ dw from the Z decay 
rate, from the forward-backward asymmetry and from LR polarization measurements, the 
top quark mass has been determined up to a residual uncertainty of less than 10 GeV 
plus an additional uncertainty of about 20 GeV induced through the variation of vtin 

mi = 178 ± 8^^^ GeV. (1.19) 

The close agreement between direct and indirect top mass determination can be con- 
sidered a triumph of the Standard Model. Its predictions are not only valid in Born 
approximation, as expected for any effective theory, also quantum corrections play an 
important role, and are indirectly confirmed. Encouraged by this success and in view of 
the improved accuracy of theory and experiment it is conceivable that the same strategy 
can lead to a rough determination of the Higgs mass, or, at least, to a phenomenologically 
relevant upper limit. 



1.1.3 The quadratic top mass dependence of bp 

The quadratic top mass dependence of bp is a cornerstone of the present precision mea- 
surements ||2^. In view of its importance and the pedagogical character of these lectures 
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it is perhaps worthwhile to present a fairly pedestrian derivation of this result. 

Let us first consider the definition of the weak mixing angle in Born approximation. It 
can be fixed through the relative strength of charged vs. neutral current couplings: 




2^2 



9 



7m(1 - 75) 



;i.2o) 



4 cos 9w 



(2// - 4g; sin^ Ow) 7m + 2//7m75] (1-21) 



with the SU(2) coupling related to the electromagnetic coupling through g = e/ sin 9w- 
Alternatively sin^ 9w is defined through the mass ratio 



sin^ 9 



w 



These two definitions coincide in Born approximation 

9w = 9w 



Born 



(1.22) 



;i.23) 



However, the self energy diagrams depicted in Fig. |L4| lead to marked differences between 
the two options, in particular if rrit ^ Mw,z ^ f^b- 




Figure 1.4: Self energy diagrams which induce the mass splitting between W and W 
and a p parameter different from one. 



This difference can be traced to a difference in the mass shift for the W and the Z boson. 
For a simplified discussion consider, in a first step, sin^ 9w = and hence the SU (2) part 
of the theory only. The neutral boson will be denoted by W^. In the lowest order this 
implies 

M_ = Mo = M (1.24) 
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and the couplings simplify to 




7m(1 - 75) 



2^2 



|2//7^(l-75) 



(1.25) 



(1.26) 



In order g"^ the propagators of charged and neutral bosons are modified by the self energies 

1 1 



1 



M2 - s-n+(s) 
1 



M2 - s M^-s- Uo{s) 

The mass shifts individually are given by 



[1.27) 



(1.28) 



They are most easily calculated through dispersion relations from their respective imagi- 
nary parts. These can be interpreted as the "decay rate" of a fictitious virtual boson of 
mass k: 



kT{WI -^t + b)^Im n+(fi:2) 
The decay rates of the virtual bosons are easily calculated (m = m^p) 



;i.29) 



Kr{wi th) 



127rlv2V2j y 2k^) 



m 



Q{K^-m^). (1.30) 



The factor 3 originates from color, the factor 2 from the identical vector and axial contri- 
butions, the squared matrix element and the phase space are responsible for the second 
and third factors in brackets respectively. 

Similarly one finds 



K 



r(W* ^ hh) + nV{W; ^tt)^ 



Utt V4 



m 



4m 



2\ 1/2 



^2 -M2 



;i.3i) 
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With the large behaviour of Imll given by Imll ~ the dispersive integral eq. ( |1.29| ) 
is evidently quadratically divergent. In the limit of large k the leading (~ k^) and next- 
to-leading const) terms of eqs. ( 1.30|) and ( |1.31|) coincide. The leading and next to 



leading divergences can therefore be absorbed in a SU{2) invariant mass renormalization. 
The relative mass shift, however, the only quantity accessible to experiment, remains 
finite and is given by 



6M^ = 5Ml-5Ml 

,2 



3 g^l dn^ . 



127r 8 TT 7o k2 - M2 

-2(l-^)(l-^-i^) xe(«:2 
X <! +1 xe(/s:2) (1.32) 

xe(«:2-4m2) 

We are only interested in the leading top mass dependence: m? ^ M^. The leading term 
is obtained by simply setting ^ in the integrand. Introducing a cutoff ^ 
the leading contributions to the three integrals are given by 

3 

-2A^ + 3m^ In — + -m^ + . . . 

2 

+A2 

-?,m'\n^\ (1.33) 
J 

and hence 

Up to the proportionality constant this result could have been guessed on dimensional 
grounds from the very beginning. 

It has become customary to express the SU{2) coupling in terms of Gp and the W mass 

f - (1.35) 



such that 

6M^ 3 



:V2GFm^. (1.36) 



M2 167r' 

The ratio of neutral versus charged current induced amplitude at small momentum trans- 
fers is thus corrected by a factor 

with Sp given in eq. ( |1.36|) . 
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To discuss the phenomenological implications of this result it is now necessary to rein- 
troduce the weak mixing between the neutral SU{2) and U{1) gauge bosons. The gauge 
boson masses are induced by the squared covariant derivative acting on the Higgs field 

+ (1.38) 
giving rise to the following mass terms in the Lagrangian 



2 



+'^g\wl + wi){i + 5p) (1.39) 



The last term has been added to represent a contribution from a non vanishing 5p, induced 
e.g. by the large top mass. The finite mass shift has been without loss of generality entirely 
attributed to the charged W boson. 

The mass eigenstates are easily identified from eq. ( |1.39|) 

Z = {cW^-sB) 

A = (cB + sWs) (1.40) 

with the weak mixing angle 6w 

9 



c = cos 9w = 



smew = ^^==^ (1.41) 



defined through the couplings. This definition is, of course, very convenient for mea- 
surements at LEP, where couplings are determined most precisely. The couplings of the 
photon and the Z boson are thus also given in terms of s. The masses are read off from 
eq. (ra) 

m\ = 

= -jig +g ) 

m'w = ^/(1 + M (1.42) 



and 



2 

''^W _ -2, 
9 ^ I 

TTLv 



1 + 5p) (1.43) 



which constitutes the standard definition of the p parameter. Alternatively one may define 
the mixing angle directly through the mass ratio 

= sin^ ew = l- Ml, /Ml (1.44) 
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The two definitions coincide in the Born approximation; they differ, however, for 6p 7^ 0: 



7t2 -2 



;i.45) 



It is, of course, a matter of convention and convenience, which of the two definitions (or 
their variants) are adopted. The choice of input parameters and observables will affect 
the sensitivity towards 5p — and hence towards m^. The observables which are measured 
with the highest precision at present and in the forseable future are the fine-structure 
constant a, the muon life time which provides a value for Gp and the Z boson mass. To 
obtain the dependence of sin^ 6^ on 5p we predict from these observables. We start 
from 



9' 1 



8 72 



and express g'^ through Gp and M| 



^2 = aV2GfMIc^ 



:i.46) 



1.47) 



Alternatively g can be related to the fine structure constant 

2 2 /-2 

g =e /s 

Note the appearance of in eq. ( [L.44D and of in eq. ( |1.41|) 
One thus arrives at 



(1.48) 



e^ = AV2GlMl^{s^ + ^5p) 



;i.49) 



or, equivalently, at 



Ana 



1 



A^/2GfMI 1 + cot2 9w5p 



2 2 
s c . 



(1.50) 



Solving for sin^ 6w (defined through the mass ratio) one obtains on one hand 



sin^ e 



w 



Ana 



\ V2GfM. 



(? — 



5p 



;i.5i) 



where the Born values for and can be taken in the correction term. The definition 
of ?,ir?6w through the relative strength of SU{2) and U{1) couplings leads on the other 
hand to 



sin^ e 



w 



Ana 



\ y^GpMl 



(? — 



bp 



;i.52) 



For the actual evaluation the running coupling a{Mz) ^ ~ 129 must be employed 
Eq. ( p..51|) and eq. (|1.52|) exhibit rather different sensitivity towards a variation of bp and 
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hence of m^, with a ratio between the two coefficients of c^/s^ ~ 3.3. For a precise com- 
parison between theory and experiment subleading one-loop corrections must be included, 
and subtle differences between variants of 9w must be taken into consideration, with 
and as most frequently used options 

With increasing experimental accuracy numerous improvements must be and have been 
incorporated into the theoretical predictions. 



The full one loop corrections including all subleading terms are know since long and 
are certainly indispensible (see 0, |^, |2^, and references therein) . 

Two loop purely weak corrections increase proportional (Gi^m^/lGyr^)^. A detailed 
discussion can be found in llTll. 



QCD corrections are available at the two- and even three-loop level p|, p3, p5 



With ct, Gpi and Mz fixed one may determine rrit either from M^r or alternatively 
from sin^ Ow (corresponding to a measurement of the left right asymmetry with polarised 
beams, the r polarisation or forward backward asymmetries of unpolarised beams). These 
measurements are in beautiful agreement with the determination of nit in production ex- 
periments at the TEVATRON (Fig. [T^p . 



Figure 1.5: 
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1.2 Top Decays 



Various aspects of top decays have been scrutinized in the hterature. The large top decay 
rate predicted in the SM governs top quark physics. Radiative correctons from QCD 
and electroweak interactions have been calculated for the decay rate and for differential 
distributions of the decay products. Non-standard top decays arc predicted in SUSY 
extensions of the SM, with t — > Hb and t — > ^7 as most promising and characteristic 
signatures. Born predictions and radiative corrections (at least in part) have been worked 
out also for these decay modes. Beyond that a number of even more exotic decay modes, 
in particular FCNC decays, have been suggested. 



1.2.1 Qualitative aspects — Born approximation 

The decay of the top quark into b + W is governed by the following amplitude 

'-t 



V2 



;i.53) 



Adopting the high energy limit (m^ > M^) for the polarisation vector e^, of the longitu- 
dinal W (corresponding to hehcity — 0) 
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the amplitude is dominated by contribution from longitudinal VF's 

1-75 



Ml = 



^b 



-t 



ig mt 71+75^ 
b — - — t 



= iV2—b{l + -f5)t 

V 



This part is thus proportional to the Yukawa coupling 
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with a rate growing proportional mf. In contrast, the amplitude for the decay into 
transverse 1^'s, is obtained with the polarisation vectors 
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;i.57) 



and remains constant in the high mass limit. The rate is governed by the gauge coupling 
g and increases only linearly with m^. The longitudinal or transversal W is produced in 
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Figure 1.6: Top decays: angular momentum conservation 



conjunction with a lefthanded b quark. The production of W^s with hehcity 
thus forbidden by angular momentum conservation (see fig. p.. 61) . 

In total one finds 



+1 is 



(/i^ = -1) : {h"^ = 0) : {h"^ = +1) = 1 : 



mr 



2M^ 



0. 



;i.58) 



The implications for the angular distributions of the decay products will be discussed 
below. The decay rate 
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increases with the third power of the quark mass and, for a realistic top mass around 180 
GeV amounts to more than 1.5 GeV, exceeding significantly all hadronic scales. Before 
we discuss the implications of this fact let us briefly pursue the close similarity between 
the coupling of the longitudinal W to the tb system and the decay into a charged Higgs 
boson in a two Higgs doublet model (THDM). The decay rate is given by (see also section 



1.2.4) 



T{t^ H+b) 



8\/27r 



1 - 



m 
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mt 



^) tan2/5 + cot2/3 
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;i.60) 



The similarity between this rate and the rate for the decay into longitudinal VT's is 
manifest from the cubic top mass dependence. The minimal value of the term in brackets 
is assumed for tan /? = 



'mt/nih. Adopting mb(running) ^ 3 GeV, the minimal value of 
the last factor amounts to about 1/30. On the other hand, in any plausible THDM the 
value of tan (3 should not exceed mt/rrib. The W decay mode will hence never be swamped 
by the Higgs channel. (This fact is of course also implied by the actual observation of 
the top quark at the TEVATRON.) Up to this point we have, tacitely, assumed the CKM 
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matrix element Vtb to be close to one. In fact, in the three generation SM one predicts 
(90% CL) 

Vtb = 0.9990 ± 0.0004 ^ BR{b) ^ 1 
Vts = 0.044 ± 0.010 ^ BR{s) ^ 0.2% 

l^td = 0.011 ±0.009 BR{d)^ 0.01% (1.61) 

on the basis of CKM unitarity. In a four generation model, however, sizeable mixing 
between third and fourth generation could arise. Methods to determine the strength of Vtb 
either through single top production at a hadron collider or through a direct measurement 
of in e~^e~ colliders will be discussed below in sects. 1273 and 3.2. 



The large top decay rate provides a cutoff for the long distance QCD dynamics. The 
implications can be summarised in the statement: "t quarks are produced and decay like 
free quarks" [^]. In particular the angular distributions of their decay products follow 
the spin 1/2 predictions. This is in marked contrast to the situation for b quarks, with 
B mesons decaying isotropically. The arguments for this claim are either based on a 
comparison of energy scales, or, alternatively, on a comparison of the relevant time scales. 

Let us start with the first of these two equivalent viewpoints: The mass difference between 
B** and B mesons amounts to 450 MeV. In the nonrelativistic quark model the B** 
is interpreted as orbitally excited bq state. With increasing mass of the heavy quark 
this splitting remains approximately constant: it is essentially governed by light quark 
dynamics. The hyperfine splitting between B* and B, in contrast, is proportional to 
the color magnetic moment and hence decreases ~ l/mq. Given a decay rate of about 
1.5 GeV it is clear that T-, T*-, and T**- mesons merge and act coherently, rendering 
any distinction between individual mesons meaningless. In fact even individual toponium 
states cease to exist. From the perturbative QCD potential an energy difference between 
IS" and 25" states around 1.2-1.5 GeV is predicted. This has to be contrasted with the 
toponium decay rate Ttt ~ ^Tt ~ 3 GeV. All resonances merge and result in an excitation 
curve which will be discussed in chapter 3. 

A similar line of reasoning is based on the comparison between different characteristic 
time scales: The formation time of a hadron from a locally produced t quark is governed 
by its size which is significantly larger than its lifetime 

^Formation ~ size ^ l/0.5GeV > mecay ~ ^/^t (1-62) 

Top quarks decay before they have time to communicate hadronically with light quarks 
and dilute their spin orientation. For sufficiently rapid top quark decay even ti bound 
states cease to exist. The classical time of revolution Trev for a Coulombic bound state is 
given by (oes = 



Trev = T^aeS^ ( 1 . 63) 

With E = —aljj-mj.f.d/2 for the ground state 

T,ev = 7r/\E\ > 1/Ttt (1.64) 
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is obtained. The lifetime of the tt system is too small to allow for the proper definition 
of a bound state with sharp binding energy. 



1.2.2 Radiative corrections to the rate 

Perturbative corrections to the lowest-order result affect the total decay rate as well as 
differential distributions. Their inclusion is a necessary prerequisite for any analysis that 
attempts a precision analysis of top decays. Both QCD and electroweak corrections are 
well under control and will be discussed in the following. 



QCD corrections 

The correction to the decay rate is usually written in the form 



Born 



1 + 5, 



QCD) 



Born 
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2a, 

3 TT 



/ 
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The correction function / has been calculated in for nonvanishing and vanishing b 
mass. In the limit ml/mf — > the result simplifies considerably, but remains a valid 
approximation (Fig.p 



/ = -^i/-^o 

J-Q = 2(l-y)2(l + 2y) 

^1 = J'o[^^ + 2Li2{y)-2Lt2{l-y) 

+Ay{l -y- 2y^) \ny + 2(1 - y)2(5 + %) ln(l - y) 
-(l-y)(5 + 9i/-6i/2), 



;i.66) 



where y = m^r/mf. In the limit mf ^ / is well approximated by f{y) = 2/37r^ 
5/2 — 3y + 9y'^/2 — 3y^ \ny ^ 4. For rrit ^ 180 GeV the QCD correction amounts to 



Sqcd ~ -3.7as/TC 



1.67) 



and lowers the decay rate by about 10%. This has a non-negligible impact on the height 
and width of a toponium resonance or its remnant. 

The corrections are presently unknown, and the scale fi in as{^^) is uncertain. Indica- 
tions for a surprisingly large correction of order a^, corresponding to a rather small scale 
have been obtained recently. Diagrams with light fermion insertions into the gluon prop- 
agator have been calculated numerically ||28| and analytically p9| in the limit rrit ^ 
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Figure 1.7: QCD correction function / for the top quark decay rate. Solid line: exact 
form; dashes: rrib = 0, dotted: approximate form for large rrit/Mw- 

The BLM prescription [^] suggests that the dominant coefficients can be estimated 
through the replacement 




(1.69) 



and absorbed through a change in the scale. For the problem at hand this corresponds to 
a scale fi = 0.122mt resulting in a fairly large effective value of as of 0.15 instead of 0.11 
for /i = rrit. 



Electroweak corrections 



Electroweak corrections to the top quark decay rate can be found in ^ . They involve 
a large number of diagrams. For asymptotically large top masses the Higgs exchange 
diagram provides the dominant contribution. Defining the Born term by means of the 
Fermi coupling Gp, one derives in this limit 



r 



r(G'F)Born 
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;i.7o) 



17 , mi, 
4 mf 



+ sub leading terms. 



While the Higgs-top coupling is the origin of the strong quadratic dependence on the 
top mass, the Higgs itself is logarithmically screened in this limit. However, the detailed 
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analysis reveals that the subleading terms are as important as the leading terms so that 
finally one observes only a very weak dependence of 5ew on the top and the Higgs masses, 
Fig.pr8|. The numerical value of the corrections turns out to be small, 5ew ~ +2%. 
Electroweak corrections in the context of the two Higgs doublet model can be found in 



33[] and are of comparable magnitude. 



Figure 1.8: QCD and electroweak radiative corrections to the top decay width; adapted 
from Ref. [ 



The positive correction 5ew is nearly cancelled by the negative correction 5t of -l.l 



from the nonvanishing finite width of the W . The complete prediction taken from 
displayed in table LI for the choice ^^(/^^ 
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mr 



For /i = 0.112m( the QCD correction 



mt [GeV] 




r^aTw. [GeV] 


4^^^ [%] 


5qcd [%] 


5ew [%] 


r [GeV] 


170 


.108 


1.41 


-1.52 


-8.34 


1.67 


1.29 


180 


.107 


1.71 


-1.45 


-8.35 


1.70 


1.57 


190 


.106 


2.06 


-1.39 


-8.36 


1.73 


1.89 



Table 1.1: Top width as a function of top mass and the comparison of the different 
approximations. 

amounts to -11.6 % instead of -8.3%. This variation characterises the present theoretical 
uncertainty, which could be removed by a full al calculation only. Additional uncertain- 
ties, e.g. from the input value of as (~ 1%) or from the fundamental uncertainty in 
the relation between the pole mass irtt and the experimentally measured excitation curve 
(assuming perfect data) of perhaps 0.5 GeV can be neglected in the forseable future. 
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Hence, it appears that the top quark width (and similarly the spectra to be discussed 
below) are well under theoretical control, including QCD and electroweak corrections. 
The remaining uncertainties are clearly smaller than the experimental error in Ft, which 
will amount to 5-10% even at a linear collider [051. 



1.2.3 Decay spectra and angular distributions 

Born predictions 

Arising from a two body decay, the energy of the W and of the hadronic system (= b jet) 
are fixed to 



El 



w 



+ - ml 
2nit 

ml + ml - mly 



Eh = ' ' ^ (1.71) 

2mt 

as long as gluon radiation is ignored. The smearing of this 6 spike by the combined effects 
of perturbative QCD and from the finite width of the W will be treated below. 

Top quarks will in general be polarized through their electroweak production mechanism. 
For unpolarized beams and close to threshold their polarization is given by the right/left 
asymmetry which would be measured with longitudinally polarized beams ||36| : 



Pt = Arl (1.72) 

For fully longitudinally polarized electron (and unpolarized positron) beams the spin of 
both t and t is aligned with the spin of the e~ . Quark polarization then leads to angular 
distributions of the decay products which allow for various tests of the chirality of the 
tbW vertex. 

The angular distribution of the longitudinal and transverse W's is analogous to those of 
p mesons from r decay (m^ mt, mp — > my/) 

^^(l^P.co..) for .„-.|(Y) (1.73) 



and, after summation over the W polarizations 

^'^-P^cos^l (1.74) 



dN _ 1 A ml - 2ml 
d cos 6 2 I ml + 2m^ 



The angle between top quark spin and direction of the W is denoted by 9. In the limit 
of mt ^ Mw the coefficient of the Pt cos 9 term rises to 1, for mt = 180 GeV, however, it 
amounts to 0.43 only. 

The angular distribution of leptons from the chain t b + i'^u) will in general 

follow a complicated pattern with an energy dependent angular distribution 

dN 

fix) +g{x)Pt cos 9. (1.75) 



dxd cos 9 
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In the SM, however, a remarkable simphfication arises. Energy and angular distribution 

dN 



factorize 36, 37 



1.76) 



dxd cos 9 

This factorisation holds true for arbitrary mt and even including the effect of the nonva- 



nishing 6-quark mass |34 



QCD corrections 



The 5 spike in the energy distribution of the hadrons from the decay t ^ h + W is smeared 
by quark fragmentation (not treated in this context). 

Hard gluon radiation leads to a slight shift and distortion of the energy spectra with a 
tail extending from the lower limit given by two-body kinematics upwards to nit — 



mi + mi 
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2mt 



< Ehad <mt- mw 



;i.77) 



Including finite W-width effects and m^ ^ the differential hadron energy distribution has 
been calculated in |^8|. The hadron energy distribution is shown in Fig. p..9| for = 180 



GeV. 

The lepton spectrum (as well as the neutrino spectrum) receives its main correction close 
to the end points where the counting rates are fairly low. 

Including QCD corrections ^ |4^, ^ the spectrum of both charged lepton and 
neutrino can be cast into the form 



dN 



dxd cos 9 



A{x) + B{x) cos 9 



;i.78) 



The shape of the charged lepton spectrum is hardly different from the lowest order re- 
sult with main corrections towards the end point. Be{x) ~ Ae{x) remains valid to 
extremely high precision |^. The charged lepton direction is thus a perfect analyser 
of the top spin, even after inclusion of QCD corrections. A small admixture oi V + A 
couplings will affect spectrum and angular distributions of electron and neutrino as well. 
Assuming a. V + A admixture of relative rate k"^ = 0.1, the functions A^, and A^, are 
only marginally modified (Fig. |1.10| and p. . 1 1| ) . The angular dependence part of the neu- 
trino spectrum B^, however, is changed significantly (Fig. |1 . 1 1|) . This observation could 
provide a useful tool in the search for new couplings. 



1.2.4 Non-standard top decays 

The theoretical study of non-standard top decays is motivated by the large top quark 
mass which could allow for exciting novel decay modes, even at the Born level. A few 
illustrative, but characteristic examples will be discussed in some detail in the following 
section. 
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Figure 1.9: Distribution of the W energy for = 180 GeV without (dashed) and with 
(sohd curve) QCD corrections. 
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Figure 1.10: The coefficient functions a) A.(^{x) and b) ^i{x) defining the charged lepton 
angular-energy distribution for y = 0.25 and as{mt) = 0.11 : — - sohd hues and 
— 0.1 - dashed hues. 




0.4 0.6 0.8 1 0.4 0.6 0.8 1 



Figure 1.11: The coefficient functions a) A,^{x) and b) Bj,{x) defining the neutrino angular- 
energy distribution for y — 0.25 and as{mt) — 0.11 : — - solid lines and — 0.1 - 
dashed lines. 
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Charged Higgs decays 



Charged Higgs states appear in 2-doublet Higgs models in which out of the eight 
degrees of freedom three Goldstone bosons build up the longitudinal states of the vector 
bosons and three neutral and two charged states correspond to real physical particles. 
A strong motivation for this extended Higgs sector is provided by supersymmetry which 
requires the Standard Model Higgs sector to be doubled in order to generate masses for 
the up and down-type fermions. In the minimal version of that model the masses of 
the charged Higgs particles are predicted to be larger than the W mass, mod. radiative 
corrections, 

m{H^) > ■m(W^) [mod. rad. corr.] 

We shall adopt this specific model for the more detailed discussion in the following para- 
graphs. 

If the charged Higgs mass is lighter than the top mass, the top quark may decay into 
plus a b quark 

The coupling of the charged Higgs to the scalar (t, b) current is defined by the quark 
masses and the parameter tan j3, 



J{b,t) 
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[{nib tan P + nit cot (3) — 75 {mf, tan P — nit cot /?)] . 



;i.79) 



The parameter tan (3 is the ratio of the vacuum expectation values of the Higgs fields 
giving masses to up and down-type fermions, respectively. For the sake of consistency, 
related to grand unification, we shall assume tan /3 to be bounded by 



1 < tan /3 
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with V 



{V2Gj 



-1/2 



corresponding to the ground state of the Standard 



Model Higgs field. The width following from the coupling ( p,.79D has a form quite similar 
to the Standard Model decay mode [see e.g. |^ ], 
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The branching ratio of this novel Higgs decay mode is compared with the W decay mode 
in Fig. |1.12| a (The behaviour is qualitatively similar for mt = 180 GeV.) In the parameter 
range eq. ( |1.80D the W decay mode is dominant; the Higgs decay branching ratio is in 
general small, yet large enough to be clearly observable [Q. The Higgs branching ratio 
is minimal at tan P = ^mt/mb ~ 6 — 8. QCD corrections to the t — » Hb mode have been 
calculated in |E6| and electroweak corrections in 



The detection of this scalar decay channel is facilitated by the characteristic decay pattern 
of the charged Higgs bosons 



+ Ut and c + s 
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Figure 1.12: (a) Branching ratios for the decays hW^ and t— hH^ in two-doublet 
Higgs models 



Figure 1.12: (b) Branching ratios for charged Higgs decays to r leptons and quarks |^ 
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Since bosons couple preferentially to down-type fermions for tan/3 > 1, 



CS) 



72 47r 
72 47r 



tan P 



tan^ j3 + cot^ jS 



;i.82) 
;i.83) 



the r decay mode wins over the quark decay mode [Fig. |1.12[ b], thus providing a clear 
experimental signature. A first signal of top decays into charged Higgs particles would 
therefore be the breakdown of /i, e vs. t universality in semileptonic top decays. 

An interesting method for a determination of tan (3 is based on an analysis of the angular 
distribution of Higgs bosons in the decay of polarized top quarks 



dA^ 
dcos6' 



cot^ P{1 + cosOh) + ml tan^ P{1 — cos 9h) 



(1.84) 



an immediate consequence of the couplings given in (|1.79|) . 



Top decay to stop 

Another exciting decay mode in supersjTumetry models is the decay of the top to the 
SUSY scalar partner stop plus neutralinos, mixtures of neutral gauginos and higgsinos 
P9| , |50[] . This possibility is intimately related to the large top mass which leads to novel 
phenomena induced by the strong Yukawa interactions. These effects do not occur in 
light-quark systems but are special to the top. 

The mass matrix of the scalar SUSY partners {tLitR) to the left- and right-handed top- 



quark components {tL,tji) is built-up by the following elements |^| 



m 



Large Yukawa interactions lower the diagonal matrix elements 



-m 



with respect to the 



common squark mass value in supergravity models, and they mix the II and Ir states 
with the strength ~ rrit to form the mass eigenstates ti,t2- Unlike the five light quark 
species, these Yukawa interactions of 0{rnt) can be so large in the top sector that after 
diagonalizing the mass matrix, the smaller eigenvalue may fall below the top quark mass, 



mi < rrit 



possible 



The decay modes 



t— > t + neutralinos 



then compete with the ordinary W decay mode. Identifying the lightest SUSY particle 
with the photino 7 (the mass of which we neglect in this estimate) one finds 



r(t-^ t7) 8^2 
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This ratio is in general less than 10%. The subsequent t decays 

i bW, or W etc. 

lead to an overall softer charged lepton spectrum and, as a result of the escaping photi- 
nos, to an increase of the missing energy, the characteristic signature for SUSY induced 
phenomena. 

Depending on the SUSY parameters however, stop decays could even be more enhanced 
if the top is heavy. Decays into strongly coupled, fairly light higgsinos could thus occur 
frequently. 



FCNC decays 



Within the Standard Model, FCNC decays like t— cy are forbidden at the tree level by 
the GIM mechanism. However, they do occur in principle at the one-loop level, though 
strongly suppressed. The suppression is particularly severe for top decays since the quarks 
building up the loops, must be down-type quarks with ml setting the scale of the decay 



amplitude, TpcNC ~ aGjpmfmt. A sample of branching ratios is given below |5^ : 




BR{t- 
BR{t- 



cg) 

C7) 



10-12 



BR{t^ cZ) 
BR{t^ cH) 



10-12 
10-^ 



At this level, no Standard Model generated t decays can be observed, even given millions 
of top quarks in proton colliders. On the other hand, if these decay modes were detected, 
they would be an undisputed signal of new physics beyond the Standard Model. From 
such options we select one illustrative, though very speculative example for brutal GIM 
breaking. It is tied to the large top mass and holds out faint hopes to be observable even 
in low rate e~^e~ colliders. 

The GIM mechanism requires all L and R quark components of the same electric charge 
in different families to carry identical isospin quantum numbers, respectively. This rule is 
broken by adding quarks in LR symmetric vector representations to the "light" chiral 



representations or mirror quarks 54 : 



vector quarks: 



tR 

bR 



U 
D 



U 
D 



R 



mirror quarks: 



J L 



tR 

bR 



Ul 
Dl 



U 
D 



J R 



Low energy phenomenology requires the masses M of the new U, D quarks to be larger 
than 300 GeV. 
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Depending on the specific form of tlie mass matrix, mixing between tlie normal chiral 
states and the new states may occur at the level ~ y^'m/M, so that FCNC (t, c) couplings 

of the order ~ ^mtradM'^ can be induced. FCNC decays of top quarks, for example, 

BR{t^ cZ) ~ fraction of % 

are therefore not excluded. Such branching ratios would be at the lower edge of the range 
accessible at e+e~ colliders. 
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Chapter 2 

Top quarks at hadron colliders 



The search for new quarks and the exploration of their properties has been a most im- 
portant task at hadron colhders in the past. The recent observation of top quarks with 
a mass of around 180 GeV at the TEVATRON has demonstrated again the discovery 
power of hadron colhders in the high energy region. Several ten's of top quarks have 
been observed up to now. The significant increase of luminosity toward the end of this 
decade will sharpen the picture. The branching ratios of the dominant decay modes will 
be determined and the uncertainty in the top mass reduced significantly. For a detailed 
study of the top quark properties the high energy collider LHC will provide the required 
large number of top events [order 10"^]. 



The main production mechanisms for top quarks in proton-antiproton collisions, Fig.2.1, 



are quark- ant iquark fusion supplemented by a small admixture of gluon-gluon fusion ||55 

gg and qq tt 

Top production at the LHC is of course dominated by the second reaction. The TV-gluon 
fusion process ||56[ 

Wg^tb 

is interesting on its own. It is about a factor 0.1 - 0.2 below the dominant reaction and thus 
well accessible at the high energy pp colliders — and perhaps even at the TEVATRON. 



2.1 Lowest order predictions and qualitative features 



The dominant Born terms for the total top cross section in 
known to be of the form p5 



47raf 
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Figure 2.1: The main production mechanisms for top quarks in pp and pp colhders 
[generic diagrams]. 



with p = Ami I ^ ~ — p being the velocity of the t quarks in the tl cm frame with 

invariant energy V^. The total pp cross sections then follow by averaging the partonic 
cross sections over the qq and gg luminosities in pp (and similarly in pp) collisions, 



(J (pp— > tt) = / dT 

JAm'^ /s 



dC {gg) 



(2.2) 



The relative enhancement of the qq cross section by about a factor 3, as evident from the 
threshold behaviour 



9 s ^ 
7 vra^ 



48 s 



(2.3) 



has to be combined with the prominent qq luminosity at the TEVATRON. As shown in 

FiK.O 



5 TEVATRON 
0.1 LHC 



(2.4) 



which implies the dominance of qq annihilation, in contrast to the situation at the LHC, 
where gluon fusion is the dominant reaction. 

A number of important features can be read off from this lowest order result: 



Since the parton luminosities rise steeply with decreasing r, the production cross 



sections increase dramatically with the energy (Fig. |2.3|) 
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Figure 2.2: Parton luminosities for TEVATRON and LHC energies. 



Figure 2.3: Production cross section for tt pairs in pp and pp colliders: Tevatron (1.8 
GeV); Tevatron II (3.6 TeV); LHC (16 TeV); SSC (40 TeV). Ref.|4 . 
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Structure functions and quark- ant iquark luminosities in the region of interest for 
the TEVATRON, i.e. for ^ between 0.2 and 0.4 are fairly well known from exper- 
imental measurements at lower energies (combined with evolution equations) and 
collider studies. The predictions are therefore quite stable with respect to variations 
between different sets of phenomenologically acceptable parton distributions. The 
near tenfold increase of the energy at the LHC and the corresponding decrease of x 
and a/t by nearly a factor of ten leads to the dominance of gluon-gluon fusion and 
results in a significantly enhanced uncertainty in the production cross section. 

With the cross section proportional to and uncertainties in as{Mz) which may 
be stretched up to ±10% one might naively expect a resulting uncertainty in the 
predicted cross section. However, the increase in the parton cross section with 
increasing is, to some extent, compensated by a decrease in the parton luminos- 
ity (with increasing for the kinematical region of interest at the TEVATRON. 
This compensation mechanism has been studied in for inclusive jet production 



(Fig. p.4| ) and applies equally well for top quark production. 



At the TEVATRON the rapidity distribution is strongly dominated by central pro- 
duction, \y\ < 1, a consequence of the balance between the steeply falling proton and 
antiproton parton distributions. At the LHC, however, a rapidity plateau develops 
gradually and the distribution spans nearly four units in rapidity (Fig. |2l5| ). 



The transverse momentum distribution is relatively flat, dropping down to half its 
peak value at around pt ~ rnt/2, again a consequence of the competition between 
the increase of the phase space factor oc /5 in the parton cross section and the steeply 
decreasing parton luminosity (Fig. |2.6| ). At the LHC the distribution will be even 
flatter and pt values around 2mt are well within reach (Fig. p.7|) . This corresponds to 
CMS energies between 0.5 and 1 TeV in the parton subsystem and extremely large 
subenergies are therefore accessible. This opens the possibility to search for the 
radiation of W, Z, or Higgs bosons in this reaction. For high energies the suppression 



of the cross section through electroweak virtual corrections (cf. sect. |2.2.3| ) is, at least 
partially, compensated by the large logarithm In s/M^r ^ h- 
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Figure 2.4: The pp-initiated jet distribution at ^/s = 1.8 TeV normalized to the 
prediction from partons with as = 0.115 (i.e. MRS. 115). The data are the CDF mea- 
surements of (P a / dErpdr] averaged over the rapidity interval 0.1 < |?7| < 0.7. The curves 
are obtained from a leading-order calculation evaluated at = 0.4. The data are prelim- 
inary and only the statistical errors are shown. The systematic errors are approximately 



25% and are correlated between different Et points. (From |]57|). 
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Figure 2.5: Rapidity distribution of top quarks in gg/qq fusion at ^/s = 16 TeV 



[441 



2.2 QCD and electroweak corrections 

The observation of top quarks has been well established during the last year. One of 
the tools to study its properties, in particular its mass and its decay modes, is a precise 
experimental determination of its production cross section and subsequent decay in the 
t —>■ b + W channel. A large deficiency in the comparison between theory and experiment 
would signal the presence of new decay modes which escape the canonical experimental 
cuts; with t — > bH~^ as most prominent example. However, the early round of experiments 
had indicated even an excess of top events when compared to the theoretical prediction for 
m,t ~ 180 GeV. This observation was difficult to interpret and the original calculations 
were scrutinized again by various authors. In particular, the resummation of leading 
logarithms and the influence of the Coulomb threshold enhancement was investigated — 
in the end, however, the prediction remained fairly stable. 

In these lectures we will, therefore, in a first step, present the results from a complete 
NLO calculation (sect. |2.2.1|) . This is supplemented by a qualitative discussion of the 
resummation of higher order leading logarithmic terms. The infiuence of the Coulomb en- 
hancement is studied in section p.2.2| , electroweak corrections are presented in sect. |2.2.3| . 



Radiation of gluons may have a sizeable effect on the aparent mass of top quarks as ob- 
served in the experiment (sect. |2.2.4 ), with distinct differences between initial and final 
state radiation. 
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Figure 2.6: The differential cross section for p + p t{t) + X witli rrit = 120 GeV/c^ and 



IJ-R = Q = ynT't + Pt at ^/s = 1.8 TeV. Tfie cross section is shown at different values of 
rapidity for (1) dashed lines: lowest order contribution scaled by an arbitrary factor (2) 
solid lines: full order calculation. (From 



2.2.1 Next to leading order (NLO) corrections and resummation 
of large logarithms 



Higher-order QCD corrections |59, 58, 31] include loop corrections to the Born terms 



and 2^3 contributions like gg — > ttg, qq ttg etc. For qq annihilation a few characteris- 



tic diagrams are displayed in Fig. |2.8| . Real and virtual initial state radiation (Fig. |2.8| a,b) 
dominate, final state radiation from the slow top quarks (Fig. |2.8| c) is unimportant, vir- 
tual gluon exchange at the ti vertex (Fig. |2.8| d) leads to the Coulomb enhancement and 



will be discussed in sect. p.2.2| . The separation between qq annihilation and qg reactions 
(Fig. depends on the choice of the so-called factorisation scale /x^ which effectively 
enters the definition of the structure functions. 

The differential as well as the total production cross section can be cast into the following 
form 

dcr = J dxidx2F[{xi, ii^)F^{x2, H^)daij (s, . . . , /x^) (2.5) 

i,j=partons 

The renormalization scale and the factorisation scale are in general identified, 
IJ'R = fJ'F ^ fJ', Si. matter of convention and convenience more than a matter of necessity. 



41 



Figure 2.7: Transverse momentum distribution of top quarks in gg/qq at y/s = 16 TeV. 
(From g.) 



The parton distributions are extracted from structure functions as measured in deep 
inelastic scattering, and the analysis has to be taylored to the order of the calculation, i.e. 
to the NLO in the present case. The integrated expressions for the total cross sections 
can still be cast into a simple form 



iij 



(P) In 



(2.6) 



where s = X1X2S and the dominant lowest-order contributions fj;j^ (p) are given by the 
parton cross sections above; in addition /^°) = fj^q = 0. The subleading higher-order 
expressions for f^p '^"^ -^^^^ 



and are given in Refs.|^, ||6^. The heavy quarks are treated 
within the on-shell renormalization scheme with rrit being the "physical" mass of the 
top quark. Outside the heavy quark sector, the MS scheme has been employed. These 
higher-order terms have to be used in conjunction with the running coupling as (/i^) and 
the gluon/light-quark parton densities evolved in 2-loop evolution equations, /i is the 
renormalization scale, identified here also with the factorization scale; typical scales that 
have been chosen are /i 



rrit and 



More technical details are discussed in 



Ref. 441. 



The lowest- and higher-order predictions are compared with each other in Fig. 2.3. In [44 
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(b) 





(c) (d) 
Figure 2.8: Initial and finite state radiation in the reaction qq 

+ . . . 



tt 





Figure 2.9: Reaction qg ttq 



it has been argued that the subdominant 2— 3 contributions add up to less than 10% 
of the dominant lowest-order results. The theoretical uncertainties of the predictions for 



the LHC due to different parton distributions were estimated about ±10% plus a 
±10% variation due to the scale ambiguity /i. The impact of the additional shift from 
the resummation of large logs arising in higher orders will be discussed below. Note that 
the "K factor", defined formally by the higher order corrections to the LO parton cross 
section, but the parton distributions and kept fixed, amounts to an ^ 50% correction 
of the Born terms. 

It is also instructive to study separate, physically distinct components of the results 
61[] . The initial state bremsstrahlung (ISGB) processes, illustrated for the gluon initiated 
reactions in Fig. |2.10| , dominate around threshold [y/T > 2mt or pt < rrit)], the case of 
relevance at the TEVATRON. The gluon splitting (GS) and the flavor excitation (FE) 
contributions become increasingly important for 3> 2mt, the situation anticipated for 
the LHC. 

Let us concentrate in the remainder of this section on the predictions for TEVATRON 
energies. Initial state radiation reduces the effective energy in the partonic subsystem, 
requiring larger initial parton energies to reach the threshold for top pair production. 
Considering the steeply falling parton distributions Fj{x) one might, therefore, expect 
a reduction of a through NLO contributions. However, the same effect is operative in 
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Figure 2.10: Generic QCD diagrams of the dominant higher order mechanisms. 



the very definition of Fj (Fig. through deep inelastic scattering, including NLO 





(a) (h) 
Figure 2.11: Definition of quark distributions, including NLO. 

corrections. In fact, without this compensation mechanism the result would not even be 
finite. However, the magnitude or even the sign of the correction cannot be guessed on an 
intuitive basis, and, not quite unexpected, even the precise form of f^^^ and f^^^ depends 
on the definition of the structure functions. The most prominent examples are the MS 
scheme where 1/e poles (plus ln(47r) — 7) from coUinear singularities are simply dropped 
[more precisely, they are combined with the corresponding singular terms which arise in 
the NLO definition of the structure function] and finite corrections have to be applied 
when comparing to deep inelastic scattering, and the DIS scheme, where Fj are defined 
through deep inelastic scattering to all orders. 

Let us illustrate the qualitative aspects in the simpler example of NLO contributions to 
the Drell Yan process. The dominance of initial state radiation in the corrections to tt 
production will allow to apply the same reasoning to the case of interest in these lectures. 
Including NLO corrections one obtains 



'^Q Xj ^ ^» 



(2.7) 



with 



9Q^s 

e'Al - z) 
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(1)/ X 

47r2 




+ 1 H — — \ 5{1 — z) + regular functions 



(2.8) 



(The quark-gluon induced reactions will not be discussed in this connection.) The plus 
prescription which regulates the singularity of the distributions at 2; = 1 arises from the 
subtraction of collinear singularities. It can be understood by considering the limit 

/ln"(l-2)\ , fln^n-z)^, , ^ln"+'e^ , , 

\ 1 — z J ^ \y I — z n + 1 J 

with the coefficient of the S function adjusted such that the integral from zero to one 
vanishes. 

Equivalently the plus-distribution can be defined through an integral with test functions 
f{z). If f{z) vanishes outside the interval (r, 1) a convenient formula which will be of use 
below reads as follows 



Jt \ I — z , Jt L — z Jo V — z 



The Born term uj'^^^ is simply given by a 5 function peak at 2; = 1, corresponding to the 
requirement that the squared energy of the partonic system r/2; x s and the squared mass 
of the muon pair r ■ s = be equal. 0{a^ corrections contribute to the 5 function 
through vertex corrections and a continuous part from initial state radiation extending 
through the range 

T<Z<\ (2.11) 

The upper limit z = \ corresponds to the kinematic endpoint without radiation, the 
requirement t < z originates from the fact that the parton luminoisities Tijirjz^ vanish 
ioYT I z > 1. The regular and the sub leading pieces of lj are process dependent, the leading 
singularity ~ (1 + z'^) (^^^^3^) is universal (and closely related to the q qg splitting 
function) and equally present in ti production. 

The suppression of final state radiation in top pair production allows to extend the anal- 
ogy to the Drell Yan process and to employ resummation techniques that were successfully 
developed and applied for muon pair production |Q. A complete treatment of this re- 



summation is outside the scope of these lectures. Nevertheless we shall try to present 
at least qualitative arguments which allow to understand the origin of these large log- 
arithms. (For a similar line of argument see [Q.) With the energies of the partonic 
reaction = C(2 - Arrit) ^ 350 - 700 GeV and the CMS energy = 1800 GeV of 
comparable magnitude it is clear that the ratio s/s will not give rise to large logs. How- 
ever, large logarithms can be traced to the interplay between the collinear singularity in 



the subprocess and the rapidly falling parton luminosity Tij (cf. eq. ( |2.7| )). This rapid 
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decrease leads to a reduction in the effective range of integration. Let us, for tlie sake of 
argument assume a range reduced from 

T <z<l (2.12) 

to 

l-5<z<l (2.13) 
and evaluate the leading term. For a constant luminosity Tirjz) j z = T one would obtain 

If the region of integration extended through the full kinematic range and with r = 
Ami/ s = 0-04 there would be no large log. For the restricted range of integration, however, 
one finds 

' j Ml-^) \ ^ = 1^2(5)^ (2.15) 

1-5 \ 1-Z 

For small 6, corresponding in practice to steeply falling luminosities one thus obtains 
large, positive (!) corrections from the interplay between J^ij^r/z) and Ugg. 

To arrive at a reliable prediction the leading terms of the form (a^ln^ . . .)" thus have to 
be included. The results are based on alternatively momentum space or impact param- 
eter techniques which were originally developed for the Drell Yan process and applied 
to top pair production in |Q. An additional complication arises from the blow up of 



the coupling constant associated with the radiation of soft gluons for s — > Amf. This 



has been interpreted in |64| as a breakdown of perturbation theory. Different regulator 



prescriptions have been advertised. In |Q a cutoff s — Amf > /ig with /Iq <^ Amf was 
introduced to exclude a small fraction of the phase space. 

The result is fairly stable for qq induced reactions with Ho/rrit chosen between 0.05 and 
0.2. The small contribution from gluon fusion, however, is sensitive towards /io which had 
to be chosen in the range between 0.2mt and O.Sm^, a consequence of the enhancement of 
radiation from gluons. A slightly different approach ( "principal value resummation" ) has 
been advocated in which circumvents the explicit /xq dependence of the result, but 



leads essentially to the same final answer (Table |2.1|) . 

The result of the improved prediction (central value) is compared to the fixed order 
calculation (with /i^ = 4m^,m^,m^/4) in fig. p.l2| . Resummation evidently increases the 
cross section slightly above the previously considered range. The history of predictions 
is shown in table p.2| , with ^/s = 1.8 TeV and rrit = 180 GeV as reference values. The 
table demonstrates that the spread of predictions through a (fairly extreme) variation of 
structure functions (DFLM vs. ELHQ) and through a variation of the renormalisation 
and factorisation are comparable — typically around ±10%. Leading log resummation 
increases the cross sections by 10 — 15%, with a sizeable sensitivity towards the cutoff 
prescription. A reduction in rrit by 5 GeV leads to an increase of a by about 0.8 pb. 
Theory and experiment, with its present result of T.GI^Iq pb and 6.3 ± 2.2 pb from CDF 
and DO respectively are thus well compatible (Fig. p.l3| ). 
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/v>-, rr^ n 
rrit [LreVJ 


1 7?; 

i /^iD 


1 on 
ioU 


/io (mm) 


4.72 


3.86 


/io (centr) 


4.95 


4.21 


/io (max) 


5.65 


4.78 


"principal value" 


5.6 


4.8 



Table 2.1: Top production cross section (in pb) for y/s =1800 GeV for different values of 
the cutoff /io (first three lines) and for the "principal value" prescription (fourth 
line). 





a [pb] 


Altarelli et al. |59| 


3.52 (DFLM) 
4.10 (ELHQ) 


Laenen et al. 


64 1 


3.5 (/i^ = 4m2) ^ 
3.8 (/i2 = m2) 
4.05 (/i^ = mV4) ^ 


> MRSD 


Resummation 




Laenen et al. 


3.86 ^ 
4.21 
4.78 ^ 


> vary /io 


Berends et al. 


m 


4.8 central value 


Berger et al. ||65| 


4.8 "principal value res. 


)7 



Table 2.2: History of predictions for the production cross section for ^/s = 1.8 TeV and 
mt = 180 GeV. 



2.2.2 Threshold behaviour 



Near the tt threshold the cross sections are affected by resonance production and Coulomb 
rescattering forces [|6^ , 
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These corrections can be estimated in a simplified 
potential picture. The driving one-gluon exchange potential is attractive if the ti is in 
color-singlet state and repulsive in a color-octet state [|6^ , 





^ tt) = 


^ crB{99- 


^ ti) 


1^(1) |2 




^ ti) = 


^ (^b{99- 


^ ti) 


|^(8)|2 




^ ti) = 


CTBiqq- 


^ ti) 


^(8) |2 



(2.16) 
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Figure 2.12: The NLO exact cross section as a function of the top quark mass for three 
choices of scale: fi = m/2 (upper sohd hne), fi = m (central solid line) and /i = 2m 
(lower solid line), and the NLO exact cross section plus the O(a^) contribution at /i = m 
(dashed line) (from []64|). 



with the correction factors (see fig. p. 141 ) given in NLO by 

The summation of the leading iras/ (3 terms to all orders results in the familiar Sommerfeld 
correction factor 

/3|^P = /3-^ (2.18) 
1 — e ^ 

For ti in the singlet configuration x = x^^^ = |^^, for octet states x = x^^^ = 

The Coulombic attraction thus leads to a sharp rise of the cross section at the threshold 
in the singlet channel, even if no resonance can be formed anymore, since the phase space 
suppression of the Born term cr^ oc /3j is neutralized by the Coulomb enhancement of 
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Figure 2.13: Comparison of experimental results for the top quark mass and production 
cross section with the theoretical prediction. 



the wave function |\E'i|^ oc as/p. In the octet channel (dominant for qq annihilation), 
by contrast, the cross sections are strongly reduced by the Coulombic repulsion which 
leads effectively to an exponential fall-off of the cross sections oc exp [— 7ras/6/3i] at 
the threshold Due to the averaging over parton luminosities the effects are less 

spectacular in pp or pp than in e^e~ collisions. 

The enhancement and suppression factors are compared to simple phase space ~ /3 in 
Fig. p.l7| . The dotted line corresponds to the phase space factor /3, the dashed line to the 
perturbative NLO calculation ( |2.17|) , the solid line to the Coulomb enhancement given 
in eq. ( |2.18| ). The predictions for the singlet, octet (= qq), and properly weighted gluon 



fusion channel are displayed in Figs. 2.17 a, b, and c, respectively. 
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n TTo, J 4/3 singlet 
2A 1-1/6 octet 



Figure 2.14: Vertex corrections from gluon exchange in the threshold region. 
2.2.3 Electroweak corrections 

Another potentially important modification which is closely tied to the Coulomb en- 
hancement originates from vertex corrections induced by light Higgs boson exchange. In 
a simplified treatment these are lumped into a Yukawa potential 

An r 

resulting in a reduction 5E of the apparent threshold, with 5E = —200 MeV for mn = 100 
GeV as characteristic example. The change in the normalisation by +10% could become 
relevant for precision measurements. The situation is quite similar to the one discussed 
for e^e~ colliders in section 3.2. 

Genuine electroweak contributions of O (aaf ) have been calculated to both the qq and 
gg^ ti subprocesses ||69|. The corrections include vertex corrections and box diagrams 



built-up by vector bosons and the Higgs boson (Fig. p.l6|) . 



Except for a small region close to the production threshold, which is dominated by the 
Yukawa potential, the corrections are always negative; they can become sizeably large, 
in particular if the top is very heavy and if the energy of the subsystem exceeds 1 TeV, 
not uncommon for tt production at the LHC In this case however, the large negative 
corrections are compensated by positive contributions from real radiation of W, Z, or H. 
The corrections for the qq and gg subprocess as functions of the parton energies are shown 
in Fig. p.l7| . The sharp increase of the corrections close to threshold for a light Higgs 
is clearly visible and, similarly, the large negative correction for large parton energies. 
After convoluting the cross sections of the subprocesses with the parton distributions, a 
reduction of the Born cross section at a level of a few percent is observed (Fig. ^.l(j| ). 

2.2.4 Gluon radiation 

Up to this point the discussion has centered around the predictions for inclusive top 
quark production. Additional ingredients for the experimental analysis are the detailed 
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Figure 2.15: Threshold enhancement and suppression factors for singlet (upper), octet 
(middle), and properly weighted gluon fusion (lower figure) for as — 0.15. Dotted line: 
Born approximation; dashed line: NLO approximation; solid line: Coulomb enhancement. 
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Figure 2.16: Feynman diagrams contributing to the electroweak radiative corrections. 
Dashed hues stand for H, Z, 7, or W. 

topological structure of the signal, the number of jets, the characteristics of the underlying 
event, and, of course, predictions for the background. This information allows to adjust 
in an optimal way experimental cuts and to measure the top quark through a kinematic 
analysis of its decay products. As a typical example the impact of gluon radiation on 
the top mass determination has been analysed recently. An idealized study has been 
performed e.g. in |]7D[. Radiated gluons are merged with the b jet from top decay or with 
the quark jet from W ^ qq if they are found within a cone of opening angle 



i? = ^(Ar7)2 + (A0)2 (2.20) 

with respect to b, q or q and if their rapidity is below \ri\ < 2. In this case the gluon is 
considered as top decay product and hence contributes to its invariant mass. If the gluon 
jet falls outside the cuts, it is assigned the rest of the event. 

Gluon radiation associated with top quark production, if erroneously associated with 
top decay, will thus increase the apparent m^. Radiation from top decay, if outside the 
forementioned cuts, will however, decrease the measured mass of the quark. The interplay 



between the two compensating effects is displayed in Fig. p.20| . For a realistic Rcut ~ 0.6 
a reduction Am of around 2 GeV is predicted. 



2.3 Single top production 

Virtual W bosons, originating from u dW splitting, can merge with bottom quarks 
from gluon splitting g ^ bb to produce single top quarks in association with fairly coUinear 
d- and b- jets (Fig. p.21| ). The interaction radius in the QCD gg fusion process shrinks 
with rising energy so that the cross section a{gg — > tt ) ~ ct^/s [mod. log's] vanishes 
asymptotically. By contrast, the interaction radius in the weak fusion process ISH] is 



set by the Compton wave length of the W boson and therefore asymptotically non-zero, 
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Figure 2.17: Relative correction to the parton cross section for 
figure: qq — > ti, lower figure: gg — > ti) (from |69]). 



250 GeV (upper 
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Figure 2.18: Relative correction to the hadronic cross section for s = (16 TeV)^ (from 

H) 



(T G'jpm^r/2n. The subprocess has to be folded with the quark-gluon luminosities 

a{pp ^t + X)= —Cug{r/z)a{ug ^t + b + d) (2.21) 

J T Z 

plus a similar contribution from dg —>■ t+b+u. The fall-off of the total cross section a{pp — »• 
tb) is less steep than for the QCD fusion processes. As a result, the Wg fusion process 
would have dominanted for large top quark masses > 250 GeV at the LHC (Fig. |2.23|) . 



For rrit = 180 GeV, the case of practical interest, single top production is about a factor 
5 below the QCD reaction. Nevertheless, as shown in Fig. p.23| , about 10^ top quarks will 
be produced at the LHC by this mechanism at an integrated luminosity of / £ = 10^p6~^. 
Also at the TEVATRON this process should be accessible with the anticipated luminosity. 

A close inspection of the various contributions to the subprocess ug ^ t + b + d reveals 
immediately that the by far dominant part of the cross section is due to b exchange, with 
the b quark being near its mass shell. Since the b quark is almost collinear to the incoming 
gluon, this cross section is logarithmically enhanced ~ ln(m^/m^) over other mechanisms. 
This naturally suggests to approximate the process by the subprocess u + b d + t with 
the 6-quark distribution generated perturbatively by gluon splitting based on massless 
evolution equations. The weak cross sections can be presented in a compact form, 

w ^ Glrr^ {s - m^,? ^G^n^ 
^ ' 27r s (s + m\ - m\) 271 
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Figure 2.19: Gluon radiation from top production (upper) and decay (lower diagrams) 
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and identically the same expressions for the C-conjugate reactions. 

Top quarks are created in u + g collisions, anti-top quarks in d + g collisions where the 
absorption of a W~ transforms a b quark to a t quark. The naive expectation from 
valence quark counting for the ratio of t/t cross sections, cr(M — > t) : o"((i — > t) ~ 2 : 1 is 
corroborated by a detailed analysis; in fact, the ratio turns out to be 2.1 for top quark 
masses of about 150 GeV. 

The remaining possibilities for single top production are Compton scattering (Fig. |2.22| (a)) 

g + b^W + t (2.22) 

and the Drell Yan process (Fig. |2.22| (b)) 

u + d^t + b (2.23) 

The predicted cross sections are too small to be of practical interest. Single top quark 
production via the dominant mechanism (Fig. p.21| ) offers a unique way for a measurement 
of the CKM matrix element Vtb and thus, indirectly of the top quark life time. As discussed 



in section |1.2.1| Vtb is strongly constrained to be very close to one for three generations 
in a four generation model Vtb may be quite different from these expectations. 
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Figure 2.20: Shift in the apparent top mass as a function of R cut (from |170| ). 




Figure 2.21: Characteristic diagram for single top production. 
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Figure 2.22: Characteristic diagrams for top production via Compton scattering (a) and 
Drell Yan process (b) 

2.4 Quarkonium production 

Both charm and bottom quarks have been discovered at hadron colhders in the form of 
quarkonium resonances J/ip and T through their distinct signals in the channel 
The search for toponium at a hadron colhder is, however, entirely useless. The broad (F ~ 
2 GeV) resonances decay with an overwhelming probability through single quark decay 
and are therefore indistinguishable from open top quarks produced close to threshold. 

The situation could be different in extensions of the SM. Decays of a fourth generation b' 

b' ^t + W (2.24) 

are suppressed by small mixing angles. Alternatively, if my < rrit + mw, the b' ^ c + W 
mode would have to compete with loop induced FCNC decays — leaving ample room for 
narrow quarkonium states. Another example would be the production of weak isosinglet 
quarks which are predicted in Grand Unified Theories. The decay of these quarks would 
again be inhibited by small mixing angles. 

Of particular interest is the search for rjf,', the ^5*0 state composed of b' and b' ||7l], ^ [73|1 . 
It is produced with appreciable cross section. Its dominant decay mode 

rib' ^ H + Z (2.25) 

is enhanced by the large Yukawa coupling, governing the coupling of the heavy quark to 
the Higgs and the longitudinal Z. For large m^/ one obtains 

r(w ^ H + Z) ml, 

r(w-77) ^ ^'-''^ 



The branching ratios as functions of are displayed in Fig. |2.24| . The complete set of 



QCD corrections for leading and subleading annihilation decay modes can be found in 
|73[| . They do not alter the picture significantly. 

It should be emphasized that the decay rjy H + Z proceeds through the axial part of 
the neutral current coupling which, in turn, is proportional to the third component of the 
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Figure 2.23: Cross sections for several mechanisms of top quark production [^. Various 
parametrizations and models of the parton densities |62| have been used. 



weak isospin. Bound states of isosinglet quarks would, therefore, decay dominantly into 
two gluon jets. 

The cross section for open bb' production at the LHC (with = 300 GeV) amounts to 
about 100 pb. The fraction of the phase space where bound states can be formed, i.e. for 
relative quark velocity f3 < 0.1, covers around 10"^ of the relevant region 



Iod/3/3 
/o 



10"^ (2.27) 



and indeed one predicts a production cross section somewhat less than 1 pb from a full 
calculation. 

For a detailed calculation of the production cross section a proper treatment of the QCD 
potential is required to obtain a reliable prediction for the bound state wave function 
at the origin. The structure of the NLO corrections for the production cross section, 
in particular of the dominant terms, bears many similarities with the result for open 
production and for the Drell Yan process (eq. |2]^). For gluon fusion the partonic cross 
section is (in the MS scheme) given by |7^, |75| 



^ s 3 M3 
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Figure 2.24: Branching ratios of rjiy for the dominant decay modes as functions of the 
bound state mass M. 
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S(l-:){Nc 



(2.28) 



where 



F{z) = 9(1 -2) 



6z(l + 

+4(i + .(l-.)-2)ln(l-.)+4(i^^^ 
/^2(z3_ _3^_2)(z3_ 2 + 2)2 ln(z) 1 



:i + ^)3(i-^) 



1 - z 



^2.29) 



and z = M"^ js. Both Born term and the virtual correction are proportional to 5(1 — z\ 
the structure of the dominant term due to gluon splitting ~ ln(l — 2;)/(l — z^ is again 
universal. Quark-gluon and quark-antiquark initiated subprocesses of order a\ can be 
found in ||7^ , [75[] . It may be worth mentioning that the structure of QCD corrections to 
light Higgs production [|76| is nearly identical to eq. p.28| . From Fig. p.25| it is evident 
that riy states with masses up to 1 TeV are produced at the LHC with sizeable rates. 
The fairly clean signature of the Z + H decay mode might allow to discover these exotic 
quarkonia and the Higgs boson at the same time. 
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Figure 2.25: Production cross section for 77^/ including QCD corrections as function of the 
bound state mass M for ^ = 16 TeV and 40 TeV. 
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Chapter 3 

Top quarks in e^e~ annihilation 

A variety of reactions is conceivable for top quark production at an electron positron col- 



lider. Characteristic Feynman diagrams are shown in Fig. |3.1| . e"'"e~ annihilation through 
the virtual photon and Z (Fig.p.l|a) dominates and constitutes the reaction of interest for 
the currently envisaged energy region. 



In addition one may also consider a variety of gauge boson fusion reactions (Fig. |3.1| b- 
d) that are in close analogy to 77 fusion into hadrons at e^e~ machines of lower energy. 
Specifically these are single top production, 

e+e~ ue'tb (3.1) 

or its charge conjugate and top pair production through neutral or charged gauge boson 
fusion 

e+e~ uuti (3.2) 

The experimental observation of these reactions would allow to determine the coupling of 
top quarks to gauge bosons, in particular also to longitudinal W bosons and Z bosons, in 
the space-like region and eventually at large momentum transfers. This would constitute 
a nontrivial test of the mechanism of spontaneous symmetry breaking. 

The various cross sections increase with energy in close analogy to 77 reactions, and 
eventually even exceed e~^e~ annihilation rates. However, at energies accessible in the 
foreseeable future these reactions are completely negligible: for an integrated luminosity 
of 10^° cm~^, at Ecm = 500 GeV and for rrit = 150 GeV one expects about one e~^e~ti 
event (still dominated by 77 fusion). At that same energy the cross sections for e^uib+c.c. 
and uutt final states are still one to two orders of magnitude smaller. 

Another interesting class of reactions is e^e~ annihilation into heavy quarks in association 
with gauge or Higgs bosons: 

e+e" ^ tiZ (3.3) 
e+e" ^ tbW- (3.4) 
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Figure 3.1: Feynman diagrams for tt or tb production. 



e+e" ^ ttH (3.5) 
e+e" ^ tbH~ (3.6) 

Two amplitudes contribute to the first reaction |^8[: The ti system may be produced 
through a virtual Higgs boson which by itself was radiated from a Z (Fig. |3.2D . The 
corresponding amplitude dominates the rate and provides a direct measurement of the 
Yukawa coupling. The radiation of longitudinal Z's from the quark line in principle also 
carries information on the symmetry breaking mechanism of the theory. The transverse 
part of the ttZ coupling, i.e. the gauge part, can be measured directly through the cross 



Figure 3.2: Amplitudes relevant for e+e — ^ ttZ and for e^e ttH. 
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section or various asymmetries in e^e'^ —>■ ti. The longitudinal part, however, could only 
be isolated with tiZ final states. For an integrated luminosity of 10^° cm~^ one expects 
only about 40 events (see sect. 3.1.7) and it is therefore not clear whether these can be 
filtered from the huge background and eventually used for a detailed analysis. 

Light Higgs bosons may be produced in conjunction with ti DTOf. They are radiated either 
from the virtual Z with an amplitude that is present also for massless fermions or directly 
from heavy quarks as a consequence of the large Yukawa coupling (Fig. p.2|) . The latter 
dominates by far and may therefore be tested specifically with heavy quark final states. 
The predictions for the rate will be discussed in sect. 3.1.7. Depending on the mass of 
the Higgs and the top quark, the reaction could perhaps be detected with an integrated 
luminosity of 10^° cm~^. 

Top quark production in 77 collisions is conceivable at a "Compton collider" . It requires 
special experimental provisions for the conversion of electron beams into well-focused 
beams of energetic photons through rescattering of laser light. A detailed discussion can 
be found in PU|. 

Chapter 3 will be entirely devoted to tt production in e^e^ annihilation. Section 3.1 will 
be concerned with the energy region far above threshold — with electroweak aspects as 
well as with specific aspects of top hadronisation. The emphasis of section 3.2 will be on 
the threshold region which is governed by the interplay between bound state formation 
and the rapid top decay. 



3.1 Top production above threshold 
3.1.1 Born predictions 

From the preceding discussion it is evident that the bulk of top studies at an e^e~ collider 
will rely on quarks produced in e^e~ annihilation through the virtual 7 and Z, with a 
production cross section of the order of cxpomt- For quarks tagged at an angle the 
differential cross section in Born approximation is a binomial in cos d 

^ (1 + cos^ 'd) ajj + - s\v? 'daL ^-"^ cos 'dap (3.7) 



dcos^ 8 V / ^ 4 4 

U and L denote the contributions of unpolarized and longitudinally polarized gauge bosons 
along the d axis, and F denotes the difference between right and left polarizations. The 
total cross section is the sum of U and L, 

a = au + ai (3.8) 

the forward/backward asymmetry is given by the ratio 

A^^ = (3.9) 
A a 
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The (T* can be expressed in terms of the cross sections for the massless case in Born 
approximation, 

4 = ^(1-/5^/3^^^ 



2 

= p^a^^ (3.10) 



with 



a 



vv 47ra^(s)e^e^ 



s 



Gf{s) , , , m|(s-m|) 



-eeCqiVe + pae)vQ- 



/ 9 9 ^ \ 9 Tri%s 



— ( 2t;eae + p(t;e + aj 1 t;QaQ ; -3 (3.11) 



The fermion couphngs are given by 

Vj = 2J/ - 4e/ sin^ , a/ = 2// (3.12) 

and the possibihty of longitudinal electron polarization (p = — 1;+1;0 for righthanded; 
lefthanded; unpolarized electrons) has been included. Alternatively one may replace 
Girml by 

^^"^^ = AT • 2fl (3.13) 

V 2 sm b'vi/ cos^ t^vy 



With sin Qy^ (^0.23) interpreted as sin 6*5// [0 this formula accommodates the leading 
logarithms from the running coupling constant as well as the quadratic top mass terms 
in the threshold region. 

3.1.2 Radiative corrections 

QCD corrections to this formula are available for arbitrary vn} j s up to first order in as- 

2 V vr 

(1 + ^^i^^) (3.14) 
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The exact result for Ky^^ can be found in |Q. These QCD enhancement factors are 



well approximated by |84 



Ka 



4 
4 



TT 

w 

TT 

w 



3 + /3 
4 
19 

To ~ 



TT 

2 

22 



3 

47r 
7 



3 

47r 



127r 



a. 



251og(4piVA2) ^^-^^^ 

The next to leading order corrections to Ky were calculated only recently p5|, The 
scale in chosen above was guessed on the basis of general arguments which were 
confirmed by the forementioned complete calculations. 



2-Kas 



For small (3 these factors develop the familiar Couloumb enhancement ~ ^^f^? compen- 
sating the phase space ~ (3. This leads to a nonvanishing cross section which smoothly 
joins the resonance region. Details of this transition will be treated in section 3.2. 

To prepare this discussion, let us briefly study the limit of applicability of fixed order 
perturbation theory. The leading terms in the perturbative expansion close to threshold 
are obtained from Sommerfeld's rescattering formula {^x 

X ^ X Bix"^ 



4 7ra \ 
3 /3 y 



K 



Som 



X 

1+2 + 



2! 



B2X^ ^ Bsx^ ^ 



4! 



X X 

IH \ 

2 12 



X 



X 



6! 



+ — ±... (3.16) 
720 5040 ^ ^ 

with Bi being the Bernoulli numbers. At first glance one might require x < 1 for the 
perturbative expansion to be valid. However, significantly larger values of x are accept- 
able. The full Sommerfeld factor Ky"'^ is remarkably well approximated by the first three 
terms of the series for surprisingly large x (only 6% deviation for x = 4!). For top quarks 
this corresponds to /3 ~ 0.13 — 0.14 and hence down to about 3 GeV above the nominal 
threshold. Upon closer inspection one also observes that the formula given in eq. ( p.l5| ) (a 
result of order as) coincides numerically well with the correction factor Ky°^ (^1 — 
which incorporates rescattering and hard gluon vertex corrections. The results presented 
in these lectures are based on the Born predictions plus 0{as) corrections. 

Initial state radiation has an important influence on the magnitude of the cross section. 
a{seff) is folded with the Bonneau Martin structure function, supplemented by the sum- 
mation of large logarithms. A convenient formula for the non-singlet structure function 
in the leading logarithmic approximation has been obtained in [0, which is a natural 
extension of a formula proposed in This leads to a significant suppression by about 
a factor 



0.5-0.6 



(3.17) 



with 6W = 1 — 5 GeV in the resonance and threshold region. The correction factor 
increases rapidly with energy, but stays below 0.9 in the full range under consideration 
(Fig.0). 
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Electroweak corrections to the production cross section in the continuum have been stud- 
ied in [^. Apart from a small region close to threshold they are negative. Relative to the 
Gf parametrized Born approximation they decrease the cross section by -6.3% to -9.3%, 
if rrit is varied between 100 and 200 GeV, niH between 42 and 1000 GeV, and Ecm fixed 
at 500 GeV. QCD and electroweak corrections are thus of equal importance (Fig. p.4| ). 

Close to threshold and for relatively small Higgs boson masses a rapid increase of these 
corretions is observed (Fig. p.5|) which can be attributed to the attractive Yukawa potential 
induced by light Higgs boson exchange. Several GeV above threshold, and for itih around 
or below 100 GeV it is more appropriate to split these corrections into hard and soft 
exchange and incorporate the latter in an instantaneous Yukawa potential ||90||. 
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Figure 3.3: Cross section for tt production, including resonances, QCD corrections and 
initial state radiation in units of apoint- 



Longitudinal polarization 

It should be mentioned that linear colliders might well operate to a large extent with 
polarized (electron) beams. The cross section for this case can be derived from (|3.11|) . 
For top quarks the resulting right/left asymmetry 

Alb^ = {aL - au)/{aL + ctr) (3.18) 

is sizable (Fig.|3.6|) and amounts to about — 0.4, reducing the production cross section with 
righthanded electrons. However, selection of righthanded electron beams decreases the W 
pair cross section even stronger, thereby enhancing the top quark signal even before cuts 
are applied. Electroweak corrections to Alh in the threshold region have been calculated 
in |2[. 
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Figure 3.4: Genuine electroweak corrections to top production in e~^e annihilation. From 
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Figure 3.5: Relative size (in percent) of electroweak corrections in the threshold region 



for rrit = 200 GeV and different Higgs masses (from |^ 
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Figure 3.6: Right/left asymmetry as function of Ecm for mt = 150 GeV. 



3.1.3 Top quark fragmentation 

The experimental analysis of charm and bottom fragmentation functions has clearly 
demonstrated that heavy quark fragmentation is hard in contrast to the fragmentation 
of light quarks. This is a consequence of the inertia of heavy particles, the momentum 
of which is not altered much if a light quark q is attached to the heavy quark Q in the 
fragmentation process to form a bound state (Qq), see e.g. |Q. At the same time soft 
infrared gluon radiation is damped if the color source is heavy. 

For rrit > 100 GeV the strong fragmentation process and the weak decay mechanism are 
intimately intertwined |Q. The lifetime t^, < becomes so short that the mesonic (tq) 



and baryonic (tqq) bound states cannot be built-up anymore. Depending on the initial top 
quark energy, even remnants of the t quark jet may not form anymore Hadrons can 



be created in the string stretched between the t and the i only if the quarks are separated 
by about 1 fermi before they decay. If the flight path 7r^, is less than 1/2 fm, the length 
of the t — i string is too short to form hadrons and jets cannot develop anymore along 
the flight direction of the top quarks. For mt ~ 180 GeV top quark energies above 1 TeV 
are required to allow nonperturbative strings between t and i. "Early" nonperturbative 
production of particles from the string between t and i is thus absent for all realistic 
experimental configurations. "Late" production from the b and b jets produced in top 
decays dominates. This is illustrated in Fig. |3.7| . Early production dominates for mt = 90 
GeV, late production for mt = 120 GeV and a forteriori for the actual value around 180 
GeV. 



The perturbative radiation of soft gluons, too, is interrupted by the t quark decay |9^ . 
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Figure 3.7: Rapidity distributions of "early" particles (full) and "late" ones (dashed) for 
three different top masses: 90 GeV in a) 100 GeV in b) and 120 GeV in c). (From ||94|| .) 



70 



Figure 3.8: Sensitivity of the soft gluon distribution to the top width for back-to-backk t 
and b: dN / d cos{e g) for M = 140 GeV, 9b = 180°, and 1 = 0, 0.7, 5, 20 GeV, and oo as 
marked. W = 1000 GeV. (From Q.) 



The angular distribution (9) and the energy distribution 
approximately given by 

e^de^ 



dPr, 



4a., 



du 



3n 



e2 + ^ 



+ 



of the radiated gluons is 



(3.19) 



for a short-lived radiation source accelerated to 7 = Et/rrit. The gluons accumulate on 

the surface of a cone with half-aperture 9c ~ for a long-lived t, but ~ •j^^^J'jT/u 
if the particle decays quickly. The energy spectrum rises from zero to a maximum at 
~ 7r before falling off ~ l/cu for large u, if the width is greater than the confinement 
scale A. 

The impact of the finite width on the angular distribution of gluon radiation will be 
visible if u < TtEt/rrit. For a linear collider with cm. energy of 2 TeV gluon jets with 



energies of 10 GeV and below would be affected |96]. The radiation pattern is shown in 
Fig. |3.8| for uit = 140 GeV and ^/s = 600 GeV, with Fj tuned to different values in order 
to demonstrate the sensitivity of such a measurement. The impact of the conversion of 
gluons to hadrons has been ignored in this study. The picture is, further, complicated 
by the interference between radiation from top production and decay — a phenomen 
characteristic for unstable particles. These phenomena allow to probe the time evolution 
of hadronisation in a unique manner. Their understanding is a necessary prerequisite for 
any top mass measurement through a kinematic analysis of 6-jet-l^ final states. 
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3.1.4 Static t parameters 



Because of the large t mass, deviations from the Standard Model may manifest themselves 
in the top quark sector first. Examples in which the large mass is crucial are provided by 
multi-Higgs doublet models, models of dynamical symmetry breaking and compositeness. 
These effects can globally be described by form factors parametrizing the electroweak ti 
production current (a = 7, Z) and the weak (t, h) decay current (a = — ) 



3l ~ F^lI.Pl + F^^^.Pn + ^ [F!^,Pl + F^^P^] (3.20) 

[Pl,r project on the left and right chirality components of the wave functions.] In the 
Standard Model, F{i^ = 1 while all other F~ vanish; F^^^ = F^^ = 1 and F21 = -F^r, — 0) 
analogously for the Z current. CV invariance requires = -^2/?^ in the tt production 
current, and equal phases for and Fg"^ etc. in the decay current. The static values of 
the form factors ^2^'^ are the anomalous magnetic and electric dipole moments of the top 
quark. 

The form factors are determined experimentally by measuring the angular distribution 
of the ti decay products, e'^e"— > ti,t-^ bW~^, ff' etc. This requires the top quark 

to be treated as a free particle, the polarization of which not being affected by non- 
perturbative hadronic binding effects. This assumption is justified by the short lifetime 
of the top quark as discussed earlier. Details of the general helicity analysis can be found 
in the literature [0, |9^ . 



3.1.5 Normal polarization of the top quarks 

A non-zero component of the t polarization vector that is normal to the production 
plane can be generated only by the interference between complex helicity flip and non- 
flip amplitudes. Such relative phases can arise from CV violation but also from higher 
order loop corrections due to gluon exchange in the final state | p.O(j| , p,01| , P7[] or electroweak 



corrections involving Higgs and gauge bosons [0. The QCD induced normal polarization 



is generally less than 5%, the electroweak normal polarization is smaller still. [By contrast, 
longitudinal and transverse polarization components within the ti production plane are 
generated already at the tree level of the electroweak interactions and they are large in 
general; see ||100| ] for the discussion of details.] 

3.1.6 Angular correlations of ti decay products 

As stated in the previous chapter, top quarks are produced through the virtual photon 
and Z. In the threshold region they are polarized to a degree 

Pt = Arl ^ -0.4 (3.21) 
Assuming for the distribution of leptons from the decay of polarized top quarks 

dN 

-— = fix) + g{x)P, cos e (3.22) 

ax a cos U 
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(with g{x) = f{x) in the Standard Model, see eq. |1.75| ) the angular distribution allows 



to test for the chirality of the tb current. Implicitly it was assumed that hadronization 
does not affect the top spin degrees of freedom |102|| . This assumption can be tested 



independently through the study of correlations between t and i decay products. In the 
threshold region the spins are correlated oc (1 + |s+ ■ s"!^). This leads to the following 
correlated distribution: 

, /a n = U{^+)f4x-) + \94x^)9-{x-) cos0+_ (3.23) 
ax +ax-a cos tlj^^ 6 

where /+ = /_ and g+ = —g- . denotes the angle between and i- . After averaging 
the lepton energies, 

^ = 1 + Ih+h. cos (3.24) 

a cos 3 

Note that the coefficient of the correlation term is —h'^/3 and hence always negative 
(assuming CP conservation). Since < 1 it ranges between and —1/3. This limiting 
value is assumed in the Standard Model. A detailed discussion with illustrative examples 



is given in ||45|| . 



3.1.7 Testing the Yukawa Coupling 

With its relatively large Yukawa coupling gy = \/2mt/v ~ 1 the top quark is uniquely 
suited to test one of the basic ingredients of the Standard Model, the coupling between 
top quarks and the Higgs boson. The verification of this crucial prediction would confirm 
the mechanism for the generation of fermion masses and hence complete the measurement 
and analysis of basic couplings. Alternatively, any deviation would provide unambiguous 
proof for new physics. 

Different strategies are at hand at an e^e~ collider which are closely tied to the cms 
energy available and to the mass of the Higgs boson. For relatively light bosons a variety 
of possibilities appear to be promising: vertex corrections affect the cross section for tt 
production in the threshold region. For a collider in its early stage with an energy around 
500 GeV it may well be the only option available and will be discussed more thoroughly 
in section 3.2. For higher energies, say around 1 TeV, a promising choice is the radiation 
of Higgs bosons from tt (see Fig. p.2|) with a cross section around 1 fb (Fig. p.9|) ||79 



Alternatively one may analyse the top quark final states in conjunction with a Z boson. 
This reaction receives important contributions from the left one of the diagrams depicted 
in Fig. |3]^ if the Higgs mass happens to be relatively close to 2mt Again, for simple 



kinematical reasons, high energies are crucial for the reaction to be accessible (Fig. |3.10|) . 



3.2 Threshold behaviour 

The previous section dealt with top quark production sufficiently far above threshold for 
the reaction to be well described by the Born cross section, modified slightly by QCD 
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Figure 3.9: The cross section cr(e+e — > ttH) (from [f7^). 

and electroweak corrections. This is in contrast to the situation in the threshold region, 
where QCD plays an important role and controls the cross section. Strong forces modify 
the Born prediction. They compensate the phase space suppression and enhance the 
production rate significantly, leading to a step function like behaviour at threshold. The 
large top decay rate also plays an important role. Quarkonium resonances cease to exist 
and merge into a structureless excitation curve which joins smoothly with the continuum 
prediction above the nominal threshold. 

This sharply rising cross section allows to study top quarks in a particularly clean envi- 
ronment and with large rates. The following physics questions can be addressed: 

• The QCD potential can be scrutinized at short distances, with the non perturbative 
tail cut off by the top decay. As a result Aqcd or as could be determined accurately. 

• The top quark mass can be measured with a precision which is only limited by the 
theoretical understanding of the excitation curve, but in any case better than 500 
MeV. 

• Top quarks are strongly polarized (about 40%) even for unpolarized beams; and 
longitudinal beam polarisation will enhance this value even further. Detailed studies 
of top decays, in particular of the V — A structure of the tbW coupling are therefore 
feasible. 

• The interquark potential is — slightly — modified by the Yukawa potential induced 
by Higgs exchange. The excitation curve and the top quark momentum distribution 
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Figure 3.10: The Higgs mass (m//) dependences of the total cross sections of e+e —>■ ttZ 
for various top quark masses nit. The cm. energy a/s is set to be 500 GeV (a) and 1 TeV 
(b). (From [^.) 
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may therefore lead to an indirect measurement of the Yukawa couphng. 

• The large number of top quarks in combination with the constrained kinematics at 
threshold, could facilitate the search for new decay modes expected in extensions of 
the SM. 

With this motivation in mind the following points will be discussed: After a brief review 
of qualitative features of threshold production (section 3.2.1), the present status of our 
theoretical understanding of the total cross section will be presented in section 3.2.2. 

The momentum distribution of top quarks and their decay products offers an alternative 
and complementary route to probe the interquark potential, as shown in section 3.2.3. 
Spin effects and angular distributions are sensitive towards the small P-wave contribution 
induced by the axial part of the neutral current. The theoretical framework and the 
resulting predictions are collected in section 3.2.4. Rescattering, relativistic corrections 
and other terms of order will be touched upon in section 3.2.5. 



3.2.1 Introductory remarks 

For a qualitative understanding it is illustrative to compare the different scales which 
govern top production close to threshold. The quarks are produced at a scale comparable 
to their Compton wave length 

(iprod ~ l/m (3.25) 

Electroweak vertex corrections do not alter this behaviour significantly, since Z- or W- 
boson exchange proceeds at a distance ~ l/mz which is still short compared to scales 
characteristic for the bound state dynamics. For the QCD potential 

4 Oi 

Vqcb = -gy (3.26) 

one anticipates an effective coupling constant ~ 0.2, ii as is evaluated at the scale of 
the Bohr momentum 

ke ~ Y ~ ^OGeV (3.27) 

The resulting Bohr radius 

rBohr = 1/kB (3.28) 
is small compared to hadronic scales. The binding energy of the IS level 

Eb= (^a,y^«2G'ey (3.29) 

and, quite generally, the separation between different resonances, is smaller than the decay 
rate 

2rt SGeV (3.30) 
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whence all resonances will merge and join smoothly with the continuum. 
The coupling strength k of the Yukawa potential 

2 

Vy = --e-'"^'' with K = V2Gf^ = 0.042 (3.31) 

r Air ^ ^ 

is comparable to the QCD coupling 4/30;^ = 0.2 in magnitude. The exponential damping, 
however, with a cutoff I/ttih -C reohr and a lower limit ttih > 65 GeV, reduces the impact 
of the Yukawa potential quite drastically. (The situation may be different in multi-Higgs- 
models: the couplings could be enhanced and, even more important, the Higgs might be 
lighter!) Furthermore, the nonrelativistic treatment is no longer adequate and retardation 
effects must be taken into consideration. 

The large top quark width plays a crucial role for the threshold behaviour, which is 
best understood in the framework of (nonrelativistic) Green's function techniques. The 
production of ti from a pointlike source (actually of extension l/m^) at x' with frequency 
u) = E is characterized by the time dependent Green's function G{f,'F,t) which is a 
solution of the time dependent Schrodinger equation with a pointlike source term 

{idt - H)G{f, t) = 6{f- f^)e-'^' (3.32) 

with 

H^^ + Vir) (3.33) 

In the problem at hand m = mt/2 is the reduced mass, f ihe relative distance between 
t and t, and the width F = 2Tt. The location of the source is at the origin r' = by 
convention, and the second argument of G will be suppressed, G{f,t) = G{r,0,t). For a 
qualitative discussion one may ignore the potential and obtains for a stable quark 



G(f,t) = -—e'^^'e-'^' 
2'Kr 

cm = ;^^e-^* (3.34) 

2m 

The corresponding current is flowing in radial direction from the source 

with a constant flux through a sphere around the origin, reflecting the conservation of 
probability. 

The width F = 2Ft is introduced in the Schrodinger equation through the replacement 

H^H-iT/2 (3.36) 

and, consequently, through the substitution E —>■ E + iT/ 2 in eq. (|3.34| ). The exponential 
damping of the flux in radial direction 

2 

J ^ ^ ^i(^/2m(E+ir/2)-y/2m(E-iT/2))r /g^g-^N 
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is most easily interpreted in two limiting cases. For E '^T the decrease 



is solely driven by the large width, with a cutoff ^ = \/2mT ^2.4 GeV for realistic 
parameters. For i5 3> F, on the other hand, the current decreases like 



j ~ %e-^^l^ (3.39) 



r 



In this latter case top quarks may travel appreciable distances, up to x^ut ~ v /T . However, 
for realisitic beam energies they hardly propagate beyond the perturbative region. 



Predictions for the Coulomb potential 

The large top decay rate restricts the range of sensitivity to the short distance part of the 
potential 

V(,) = -H^^ (3.40) 

which is approximately Coulombic, with a logarithmic variation of ay- Most of the 
qualitative features of top quark threshold physics can be understood even on the basis 
of the results for constant cty, which are available in analytical form. 

A remarkable feature of heavy quark production is the sharp rise of the cross section 
at threshold, a consequence of the attractive Coulomb force. The step function joins 
smoothly with the smeared resonances. 

Let us try to quantify this aspect with the help of simple nonrelativistic quantum mechan- 
ics {a = |q;5, m = ruit/l). The narrow resonances below the nominal threshold (i?=0) 
are located at En — —ERyd/n^ with Enyd — k%/{2m) — 

The production amplitude from a pointlike source with frequency cu — E is proportional 
to 

(f=0|V'n)=^n(0), (3.41) 

the rate correspondingly to 

i?~|Vn(0)|'5(£;-K) (3.42) 

The wave function at the origin decreases with the third power of the radial quantum 
number, 

|^,(0)|2^1/n3 (3.43) 
their spacing becomes increasingly dense 

AE = En^, -En = Enyd - ^^^) - ^^.^ ^^'^^^ 
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such that their average contribution to the cross section approaches a constant value 



(Fig. |3.11| ). Above threshold one has to project the state (r = 0| onto the Coulomb wave 
functions in the continuum ipsi'f^- These replace the plane waves which are appropriate 
for the case where final state interaction is absent. 

The production amplitude is thus governed by 

{f=0\ijE) = MO) (3.45) 

and the rate ~ |?/'£;(0)p. The threshold phase space factor v is thus compensated by the 
1/v singularity in |'?/'£;(0)p, and the cross section approaches a constant value for E ^ 0. 




Figure 3.11: Schematic representation of resonances (wiht the area of the boxes adjusted 
to represent the weights of the respective delta functions) and the continuum cross section. 

The explicit calculation yields 

Qtt 3 1 / 1 \ 

RiE) = —,TAE - E„) = 3g?- E ^.^-^ [E - Eny.^) (3-46) 

for energies below threshold, and 

with X = kB/k = na/P for energies above threshold.^ 

The perturbative expansion breaks down in the limit (3^0. The first term of this formal 
series 

^For a textbook discussion of Coulomb scattering states and a derivation of this result see e.g. [ |103| , |104| 
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underestimates the exact result by a factor two. 

Eq. ( |3.46| ) allows to connect the formalism based on narrow individual resonances with a 
formulation which is tailored to the situation at hand, namely wide overlapping resonances 
which merge into a smooth continuum. 

Instead of summing the contributions from a large number of high radial excitations one 
may directly calculate the imaginary part of the Greens function for complex energy 



a[e e 



tt) 



247r^a>„(s) 



1 



mi 



3 tt/V En-E-iVt 



s 



1 



16 Us 



/mG(0,0,E + zr) 



(3.49) 



The factor Pv{s) incorporates the contributions from the intermediate photon and Z and 
is given by 



s) 



VtVeMi 



y^s-Mi + iMzTz 



(3.50) 



Vf = 2/3J — 4e/ o-f = 2/3/ y = IQ sin 9y/ cos 9w 

{a = tteff = 1/128 has been adopted in the numerical evaluation. Radiative corrections 
to this formula have been discussed in [0.) 

The problem can be solved in closed analytical form for an exact Coulomb potential ||105| 



ki,2 



mt 



4:71 



k2 2kB k2 

1 arctan — 

rrit mt ki 



+ E 



TtkBn + k2[n^^jE^ + T^t+^^ 



n=i 



E + ^ 



E^ + T^tTE)/2 



1/2 



-asmt. 



+ n 



(3.51) 



To arrive at a realistic prediction of the total (and, in sect 3.2.3 of the differential) cross 
section the Coulomb potential must be replaced by the realistic QCD potential. 



3.2.2 The QCD potential 



On the basis of earlier conceptual work in | |106| , |107|| the asymptotic behaviour of the 
static potential has been derived in ||108| , |109|| . In momentum space the potential reads 
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in the MS subtraction scheme 



167r ajTsiQ^ 

3 g2 



31 

y 



9 



(3.52) 



The renormahzation scale /i^ has been chosen as Q^, and refers to the number of 
massless quarks. Employing standard arguments based on the renormahzation group, the 
expansion of o.jjgiQ'^) is given by 

«M5(<5^) _ 1 



4vr h, log (gVAy 



^ fe, iogiog(QVAy 

^0 logfgvAy 



(3.53) 



11 nt, 

3 ^' 



102 --n. 
3 ^ 



The leading behaviour of the potential at small distances [7^; = 0.5772 

167r 



V{r) 



(33-2n/)rlogl/(A- 



MS' 



1 - 



61 loglogl/(A' 



MS' 



hi logl/(A. 



M5' 



logl/(A 



MS' 



(3.54) 



is thus directly given in terms of the QCD scale parameter A. The exploration of V{r) 
for small distances could thus lead to a direct determination of A. For quark masses 
above 50-100 GeV the ground state properties become independent of the potential in the 
nonperturbative region. As discussed in the previous section the large decay rate acts as a 
cutoff and the predictions are fairly insensitive to the actual regularisation. However, an 
additional constant which can be traded against a shift in must be carefully calibrated. 

In practice one connects the theoretically predicted short distance part smoothly with the 
empirically determined potential above ~0.1 fermi. The asymptotic form given in (|3.52| ) 
is based on the assumption that n/ species of light quarks, taken as massless, contribute 
to the vacuum polarization, and heavier ones are ignored. The value of Ajg^ in ( 3.53| ) 
must be properly related to A;^ as determined from other experiments with a different 
number of effective light flavors ||110| , |111|| . For the momentum range of around 15 GeV 
explored by the tt system Uf =5 seems adequate. 

In the subsequent discussion the Green's function will be calculated in momentum space 
with the help of the Lippman-Schwinger integral equation | 112 |. The representation of the 
QCD potential in momentum space with the large behaviour given by eqs. ^.52| and 
|3.53| will be employed. The intermediate and small momentum dependence will be based 
on Richardson's potential. This choice allows to vary a-MsiM^) (or equivalently Ajj^) 
between 0.11 and 0.13, while maintaining a smooth dependence of aviQ"^) (Fig. |3.12|) . 
An additive constant in coordinate space (corresponding to a 5-function in momentum 



space) is adjusted to fix V{r = IGeV 



-1^ 



1/4 GeV for arbitrary aj-fg- This constraint 



avoids the unmotivated and uncontroled variation of the long distance part of V{r) 
a change in aj^. The potential in coordinate space is shown in Fig. p. 13 . 



with 
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10 1 10 10^ 10^ 

q [GeV] 



Figure 3.12: for different values of as{Mz): solid: 0.12, dashed: 0.11, dashed-dotted: 
0.13, dotted: 0.10 and 0.14. 
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QCD-Potential VJKT(r) [GeV] 




10 1 10 J. [GeV**(-l)] 

Figure 3.13: QCD potential in the position space Vjkt for different values of as{Mz): 
solid: 0.12, dashed: 0.11, dashed-dotted: 0.13, dotted: 0.10 and 0.14. 



83 



3.2.3 Realistic predictions for a^' 



For a realistic QCD potential the Green's function can only be calculated with numerical 
methods. An elegant algorithm for a solution in coordinate space has been suggested in 
113| |. As a consequence of the optical theorem (see also eqn. ( |3.49| )) only the imaginary 



part of G{r = 0, = 0, E+iTt) is needed to predict the total cross section. The differential 
equation for the Green's function 




G{ry = 0,E + iTt) = S{r) (3.55) 



is solved in a way which provides direct access to ImG(r = 0, = 0, + iTf) without the 
need to calculate the full f dependence. Alternatively, in ||112| , |1 14|| the Green's function 



in momentum space was obtained from the Lippmann-Schwinger equation 
G{p,E + trt) = Go{p,E + irt)+Go{p,E + trt) 
X / j^^V{p-q)G{q,E + tT,) 

Go{p,E + zTt) = ^.^ (3.56) 

The total cross section is in this case obtained from the integral over the differential 
distribution 

^ = ^p.(.)r,|G(p, E + ^ror (3.57) 

This second formulation is particularly suited to introducing a momentum and energy 
dependent width r{p,E) which allows to incorporate the phase space suppression and 
certain rescattering corrections to be discussed below in section 3.2.6. 

It is well known that the coupling of the virtual photon to the quarkonium boundstate is 
modified by "hard" gluon exchange. The vertex correction to the vector current produces 
an additional factor (^1 — for the quarkonium decay rate into e~^e~ through the 

virtual photon or Z. This factor can be calculated by separating the gluon exchange [ 115| 



correction to the vertex into the instanteneous potential piece and a remainder which 
is attributed to gluons with high virtualities of order mt. A similar approach has been 



developed in |^0|] for Higgs exchange. The vertex correction is again decomposed into a 



part which is given by the instantaneous Yukawa potential 

Vvukir) = (3.58) 

r 

with K = \/2Gml/A'K and a remainder which is dominated by highly virtual Higgs ex- 
change. The rapid increase of the correction in the threshold region (cf. sect. 3.1.2) is 
driven by the potential; the remainder, the hard vertex correction, is fairly energy inde- 
pendent. The total cross section is thus sensitive to the top mass, the width (which in the 
SM is uniquely determined by m^), the strong coupling constant and the mass of the 
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Higgs boson. This dependence is illustrated in Figs. |3.14[ - p.l7[ Apart from the trivial shift 
of the threshold due to a change in rrit the shape of a is affected by the rapidly increasing 
width of the top quark which amounts to 0.81 GeV, 1.57 GeV and 2.24 GeV for nit = 150 



GeV, 180 GeV, and 200 GeV respectively. This is demonstrated in Fig. |3.14| : A fairly 
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Figure 3.14: Total cross section as function of E = ^/s — 2mf for three values of the top 
quark mass. 



pronounced IS peak is still visible for rrit = 150 GeV, for rrit = 200 GeV, however, only 
a smooth shoulder is predicted. The behaviour is qualitatively very similar, if we keep 
rrit fixed say at 180 GeV and decrease or increase Tt be the corresponding amount. The 
shape of the cross section will therefore allow to determine the width of the top quark. 
A qualitatively very different response is observed towards a change in (Fig. |3.15|) . 
The binding energy increases with a^, the aparent threshold is thus lowered (This is the 
reason for the strong correlation between and rrit in the experimental analysis based 
on atot only |116|] .) and the height of the "would-be resonance" is increased. Even 
several GeV above threshold one observes a slight increase of the cross section with a^, a 
consequence of the enhanced attraction between t and t (cf. eq. |3.47| ). The impact of the 



running of on the shape of the cross section is evident from Fig. 3.15. The full QCD 
prediction with running (for a;^(M|) = 0.125) is compared to the prediction for a 
Coulomb potential with fixed. It is impossible to describe the height of the peak and 
the continuum above with the same value of a^, even allowing for an arbitrary additive 
constant Vq. The infiuence of a variation in mn is shown in Fig. |3.17| . Cross section 
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Figure 3.15: Total cross section as function oi E = y/s — 2mt for three values of a^. 



measurements with a precision better than 10% will become sensitive to the effect of a 
light Higgs boson. 

Up to this point the amplitude induced by virtual Z and 7 are included in Born approxi- 
mation only. Electroweak corrections and initial state radiation are neglected. A detailed 
discussion of electroweak corrections to the cross section and the left right asymmetry in 
the context of the SM can be found in |^|. The corresponding discussion for the two- 
Higgs-doublet model is presented in ||117|| . In this model one might encounter enhanced 
Yukawa couplings which would amplify the effect under discussion. 

Initial state radiation leads to a fairly drastic distortion of the shape of the cross section, 



in particular to a smearing of any pronounced structure. This is illustrated in Fig. |3.18 
where the predictions with and without initial state radiation are compared for otherwise 
identical parameters. 

Beamstrahlung and the energy spread of the beam lead to a further smearing of the 
apparent cross section. These accelerator dependent issues are treated in more detail in 
|]35|1. (For a related discussion see also ||116||.) 
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Figure 3.16: Comparison between the predicted cross section for constant (dashed and 
dotted hnes) and running (sohd hues) a^. 
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3.2.4 Momentum distributions of top quarks 

The Green's function in momentum space and the momentum distribution of top quarks 
(and thus their decay products) are intimately related. For a narrow quarkonium reso- 
nance with orbital quantum number n the quarks' momentum distribution is evidently 
given by the wave function in momentum space 



3 



For J/ip oi T this distribution is not directly accessible to experiment since these states 
decay through QQ annihilation only. For toponium, however, which is dominated by 
single quark decay, the decay products carry the information of their parent momentum 



and hence allow for the reconstruction of the original quark momentum distribution 

For one individual resonance this leads to the differential ti cross section (without Z 
contribution and transverse gluon correction). 



Once Fj is sufficiently large, interferences between different radial excitations become 
important and the right-hand side of this equation has to be replaced by the square of 
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Figure 3.17: Total cross section as function oi E = ^/s — 2mt for different values of the 
Higgs mass. 



the Green's function [nS, 112] 



^{p, E) = \G{p, E + ^T,)f (3.61) 
dp irsmf 

with 

G{p, E + iTt) = J dfe'^G{r, = Q,E + iVt) (3.62) 

As discussed in sect. 3.2.2, the Green's function can be obtained in momentum space 
as a solution of the Lippmann-Schwinger equation. For an energy close to the IS peak 
it exhibits a fairly smooth behaviour reminiscent of the IS wave function in momentum 
space (Fig. p.l9| ). With increasing energy an oscillatory pattern of the amplitude is 
observed, and a shift towards larger momenta (Fig. p.l9[ ). These results are intentionally 
displayed for irtt = 120 GeV, where the oscillations are still clearly visible, in contrast to 
mt=180 GeV where all oscillations are smeared by the large width Ff. The corresponding 
predictions for the distributions at rrit = 180 GeV are displayed in Fig. p.20| . The 
transition from a wide distribution below the nominal threshold to a narrow one with the 
location of the peak determined by trivial kinematics is clearly visible. The impact on 



the energy distribution of the W^s from top decay is shown in Fig. |3.21 



To characterize the momentum distribution by a single parameter, one may either choose 
its peak value, or the expectation value of the modulus of the momentum (p), the latter 
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Figure 3.18: Comparison of the tt production cross section without (dashed) and with 
(sohd hne) initial state radiation. 



being well adopted to the experimental analysis. In the situation at hand the definition 
of (p) has to be introduced with some care. The free Green's function Go (see eq. |3.56| ) 



drops ~ p~ for large momenta and this behaviour is recovered also in the presence of 
interaction. The expectation value (p) diverges logarithmically with the cutoff. In the 
narrow width approximation one finds for the leading terms 



Itdpp\Go{p,E + zT)f 
Jtdp\Go{p,E + zT)f 

T\npl/{VWTT^m,] 

nitE 1 + 



E 



TT 



(3.63) 



where a cutoff Pm has been introduced. As a consequence of the small numerical prefactor 
of the divergent term and its logarithmic cutoff dependence the result is fairly insensitive 
to the exact value of the cutoff for pm of order mt. Alternatively one may replace the 
phase space element dp/mt by the relativistic version dp/E = dp I \J + p^ to obtain a 
convergent result. In future measurements the cutoff will presumably be provided by the 
experimental analysis. 

In order to study the dependence of {p) on the strong coupling constant, consider for the 
moment the predictions for a stable quark. Some intuition and qualitative understanding 
can already be gained from the predictions based on a pure Coulomb potential ||120 . 



For a stable top quark of fixed mass the "effective threshold" can be associated with the 
location of the IS" resonance y/sthr = + Eis with Eis = —Ejiyd = —a'^nit/i which 
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Figure 3.19: Real (dashed) and imaginary (dotted) parts of the Green's function for an 
energy corresponding to the IS peak (upper figure) and for = (lower figure). Solid 
curve: |pG(p)p • 0.002. 
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Figure 3.20: Momentum distribution of top quarks for three different cms energies. 



decreases with increasing a. The height of the resonance cross section is proportional to 
the square of the wave function at the origin and hence proportional to a"^, as long as 
the resonances are reasonably well separated. In the limit of large Tt, i.e. far larger than 
Ejiyd, the overlapping 15", 2S . . . resonances have to fill the gaps between the peaks. Since 
these gaps themselves increase proportional to a^, one is left in the extreme case of large 
width with a cross section linear in a. Note that this corresponds to the behaviour of the 
cross section close to but slightly above the threshold which is also proportional to a. 

For realistic top masses one thus observes a dependence of the peak cross section linear in 
a. Since the location of the peak itself depends on a, only the analysis of the full shape 
allows to extract the relevant information. 

In a next step also the momentum distribution of top quarks has to be exploited to obtain 
further information. The discussion is again particularly simple for the Coulomb potential 
V{r) — —a/r and provides a nice exercise in nonrelativistic quantum mechanics. The 
average momentum, in units of the Bohr momentum amt/2, can be written in terms of a 
function /(e) which depends only on one variable e = E/Ejiyd if the energy E = sfs — lmt 
is measured in terms of the Rydberg energy. 

(P) = ^/(^) (3.64) 

For positive arguments the function / can be derived from obvious kinematical consider- 
ations. 
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Figure 3.21: Energy distribution of Ws from top quark decay for three different cms 
energies. 



/(e) = y/e for e > (3.65) 

For the discrete negative arguments e„ = — corresponding to the locations of the 
bound states, the radial wave functions in momentum space are given in terms of the 
Gegenbauer polynomials 

with 

;Wp)r = l. (3.67) 



(27r)= 

Using the explicit forms of 

C^{z) = 1, Cl{z) = 2z, Cl{z) =42^-1 (3.68) 
one obtains through straightforward calculation 



. I, /(_l/4) . ^, /(-1/9) . ^ 



(3.69) 
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For arbitrary n one derives the general result 



, 1 \ 8n 
^ ^ ^ (2n- l)(2n + l)7r 



with the asymptotic behaviour 



(3.71) 



This is in accord with the result expected from classical mechanics: For the average 
momentum of a particle on a closed orbit in the Coulomb potential one derives 

Quantum mechanical orbits with angular momentum zero and high radial quantum num- 
bers correspond to classical motions with excentricity e = 1 (i.e. straight lines). In this 
limiting case the classical expectation value is easily evaluated, and for (3 = 1/2 one 
finds agreement with the quantum mechanical result. For small negative energies one 
therefore obtains the behaviour /(e) = 2y/^/7r. Significantly below threshold, however, 
the average momentum obtained from the Green's function increases more rapidly with 
decreasing energy and between the IS* and the 25" state one observes an approximately 
linear dependence on the energy. 

From these considerations the dependence of the average momentum on a (with E fixed) 
is easily understandable, in particular the seemingly surprising observation that well below 
threshold (p) decreases with increasing a. From ( p.64|) one derives for a shift in a (keeping 
the energy E fixed) the following shift in (p) 




(3.73) 



Above threshold as well as close to but below threshold / oc y |e|. Hence ef'/f = 1/2 and 
the average momentum remains unaffected. The location of the minimum is thus an ideal 
place to fix the mass of the top quark. Significantly below threshold, however, ef'/f ~ 1 
and the factor in front of Sa/a becomes negative. This explains the decrease of (p) with 
increasing a. 



These results are illustrated in Fig. 3.22 . In Fig. |3.22| b we demonstrate that (p) as evaluated 



with the program for the Green function (solid line) coincides perfectly well with the 
values calculated from the analytical formula on resonance, indicated by the triangles. 
The prediction from classical mechanics, namely {p) oc -y/je[ is shown by the dotted line 
and agrees nicely for positive and negative energies. In Fig. p.22| a as is increased from 0.20 



to 0.24 and (p) changes in accord with the previous discussion. 

For definiteness we have chosen m = m(/2 = 60 GeV for the reduced mass and a = 4as/3 
with as varying between 0.20 and 0.24. The curves demonstrate the decrease of (p) by 
about 10% for the corresponding increase in a. The triangles mark the locations of the 



resonances and the expectation values for the momentum as calculated from (|3.70|) . 
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Figure 3.22: a) Average momentum as a function of E for different values of as- The 
markers sliow tlie results of the analytical calculation at IS, 2S, 3S energies, b) Compar- 
ison with the analytical result for discrete energies and with the square-root dependence 
close to threshold. 
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Figure 3.23: Energy dependence of (p), the average t quark momentum for = 0.13 
(dotted) 0.12 (dashed) and 0.11 (solid) line for rrit = 120 GeV. a) = 0.03 GeV and b) 
Ft = 0.3 GeV. 
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The qualitative behaviour remains unchanged for reahstic QCD potentials corresponding 
to different values of ^^(M^). Qualitatively the same behaviour is observed as in Fig. p.22| . 



In Fig. 3.23 a the top quark width has been set to an artificially small value of 0.03 GeV, 
in Fig. |3.23| b the realistic value of 0.3 GeV has been adopted. The finite width leads to an 
additional contribution to the momentum of order y/Tm. 



An important feature is evident from Fig. |3.23| : The momentum calculated for positive 
energy is nearly independent from ag and reflects merely the kinematic behaviour, just as 
in the case of the Coulomb potential. This is characteristic for the choice of a potential 



114 1 where the large distance behaviour is fixed by phenomenology and decoupled from 



the short distance value of a^. 

The different assumptions on the long distance behaviour are reflected in differences be- 



tween the predictions of [ p.21| , |119| , |112|| for the precise location of the ti threshold for 
identical values of ag and rrit and in differences in the ag dependence of the momen- 
tum distributions for fixed rrit and energy (see also ||116|| ). All these differences can be 
attributed to the freedom in the additive constant discussed before. The same additive 
constant appears in bb spectroscopy, such that the mass difference between top and bottom 
is independent from these considerations. 

In Fig. |3.24| the predictions for (p) vs. energy are presented for the case of a realistic QCD 
potential, assuming nit = 120/150/180 GeV. The strong rise of {p) as a consequence of 
the strong increase of is clearly visible. 



3.2.5 Angular distributions and polarization 

Close to threshold the production amplitude is dominantly S-wave which leads to an 
isotropic angular distribution. The spin of top quarks is alligned with the beam direction, 
with a degree of polarization determined by the electroweak couplings, the beam polar- 
ization and the mass of the top quark, but independent of the production dynamics, in 
particular of the potential. 

Small, but nevertheless experimentally accessible corrections do arise from the small ad- 
mixture of P-wave contributions and from rescattering of the top quark decay products. 
Let us concentrate for the moment on the first mechanism. P-wave amplitudes are pro- 
portional to the top quark momentum. For stable noninteracting particles the momentum 
vanishes at threshold. However, as discussed in the previous section the expectation value 
of the quark momentum is nonzero for all energies — a consequence of the large top decay 
rate and the uncertainty principle. Technically the P-wave contribution is calculated with 
the help of the Green's function technique. The generalization of the Lippman-Schwinger 
equation (|3.56| ) from S- to P-waves reads as follows 

np,E) = Go{p,E) + Go{p,E)J-0^^^V{p-k)J^{k,E) (3.74) 

It is then straightforward to calculate the differential momentum distribution and the 
polarization of top quarks produced in electron positron annihilation. Let us recall the 
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Figure 3.24: Energy dependence of the average top quark momentum for n 
120/150/180 GeV. Triangles, stars and circles correspond to (pf) for S-states with Vt 
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following conventions for the fermion couplings 

Vf = 2Pf - 4g/sin^ 6'w, 



21}. 



(3.75) 



P± denotes the longitudinal electron/positron polarization and x = {P+~P-)I i)-~P+P-) 
can be interpreted as effective longitudinal polarization of the virtual intermediate photon 
or Z boson. The following abbreviations will be useful below: 



ai 
a2 
as 

d 



QeQt + i^e + o-e)^td + 2qeqtVeVtd 



2veaeVfd'^ + 2qeqtaeVtd 



AvetteVtatd^ + 2qeqtaeatd 

2(^^e + al)vtatd^ + 2qeqtVeatd 
1 s 



(3.76) 



16 sin^ 6'w cos^ 6'w s - M| 



The differential cross section, summed over polarizations of quarks and including S'-wave 
and S-P-interference contributions, is thus given by 



d^ 



dp^ Anmf 



P^P- 



,ai + xa2) 1 - 



16Q;e 



+ (03 + xa4) 1 



12as\ P 



3n 



—Re ( G(p, E)F* (p, E) ) cos i9 
mt 



(3.77) 



The vertex corrections from hard gluon exchange for S'-wave and P-wave amplitudes are 
included in this formula. It leads to the following forward-backward asymmetry 



with 



^fb(p,^) =Cfb(x) ^h(p.^)' 



Reip, and 



<^(p, E) 



2 ai + XO'2 
l-4as/37r) p P*(P, ^) 



;i-8as/37r) mt G*{p,E)' 
This result is still differential in the top quark momentum. Replacing if{p, E) by 



HE) 



;i - 4«,/37r) dp ^ P*(p, E)G(p, E) 



(l-8«s/37r) /P'"dpp2|G'(p,E)r 



(3.78) 
(3.79) 

(3.80) 
(3.81) 



one obtains the integrated forward-backward asymmetry again. Again, the cutoff p^ must 
be introduced to eliminate the logarithmic divergence of the integral. 

Polarization ^12^ , \lF^ 

To describe top quark polarization in the threshold region it is convenient to align the 
reference system with the beam direction (Fig. p.25| ) and to define 



He- X nt 



rip 



X 



s± = Sn X sii. 



(3.82) 
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l3" Figure 3.25: Definition of tlie spin 

directions. Tfie normal compo- 
nent Sn points out of tlie plane. 



In the hmit of small l3 the quark spin is essentially aligned with the beam direction apart 
from small corrections proportional to (3, which depend on the production angle. A system 



of reference with sy defined with respect to the top quark momentum ||124|| is convenient 
in the high energy limit but evidently becomes less convenient close to threshold. 

Including the QCD potential one obtains for the three components of the polarization 



V^{p,E,x) = C^{x)VniP,E) sin^? 



(3.83) 
(3.84) 

(3.85) 



f3 861 with 

^ . ^ 1 04 + X«3 ^ . X I a3 + xai n ( \ 

C± X = -7^ —— , C'n(x) = -7^ —- = -^fb(x), 

2 ai + x«2 2 ai + x^i 

= Im(y9, and ip{Y>iE) as defined in ( |3.8(J| ). The momentum integrated quantities are 
obtained by the replacement ip{Y>)E) ^{E). The case of non-interacting stable quarks 
is recovered by the replacement $ — > /5, an obvious consequence of (|3.81| ). 

Let us emphasize the main qualitative features of the result: 



• Top quarks in the threshold region are highly polarized. Even for unpolarized beams 
the longitudinal polarization amounts to about —0.41 and reaches ±1 for fully po- 
larized electron beams. This later feature is of purely kinematical origin and inde- 
pendent of the structure of top quark couplings. Precision studies of polarized top 
decays are therefore feasible. 

• Corrections to this idealized picture arise from the small admixture of P-waves. 
The transverse and the normal components of the polarization are of order 10%. 
The angular dependent part of the parallel polarization is even more suppressed. 
Moreover, as a consequence of the angular dependence its contribution vanishes 
upon angular integration. 

• The QCD dynamics is solely contained in the functions or $ which is the same for 
the angular distribution and the various components of the polarization. (However, 
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this "universality" is affected by the rescattering corrections.) These functions which 
evidently depend on QCD dynamics can thus be studied in a variety of ways. 



The relative importance of P-waves increases with energy, $ ~ ^^E jmt. This is 
expected from the close analogy between $r = Re$ and (5. In fact, the order 
of magnitude of the various components of the polarization above, but close to 
threshold, can be estimated by replacing $r vl°^t- 



The Ci are displayed in Fig. |3.26| as functions of the variable x (sin 6'w = 0.2317, mt = 180 
GeV). As discussed before, C^^ assumes its maximal value ±1 for x = and the coefficient 
C|| is small throughout. The coefficient C± varies between +0.7 and —0.5 whereas Cn is 
typically around —0.5. The dynamical factors $ are around 0.1 or larger, such that the 
P-wave induced effects should be observable experimentally. 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 
-1 







: \ po 












- , , , , 1 , , , 





-1 -0.5 



0.5 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 
-1 



-1 



-0.5 



0.5 



Figure 3.26: The coefficients (|3.86|) for ^/s/2 



180 GeV. 



The normal component of the polarization which is proportional to (f-^ has been predicted 
for stable quarks in the framework of perturbative QCD [ |125| , |124|| . In the threshold 
region the phase can be traced to the ti rescattering by the QCD potential. For stable 
quarks, assuming a pure Coulomb potential V = —Aas/3r, the nonrelativistic problem 
can be solved analytically ||126|| and one finds 

hm (E-^ + zT)jG{p,E) = exp (^) T{1 + zk/p) (3.87) 
hrn (E~^ + zT)jF{p,E) = (l-t-\exp(^)T{l + tk/p), (3.88) 

£;-,p2/mt \"^* / \P/\P/ 

with k = 2mtas/3 and hence 

2 1 - AajSTT 

fAP^E) ^ -a,- — — 3.89 

3 1 — Sas/Sn 
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Figure 3.27: Real (solid) and imaginary (dashed) part of the function y9(p, E) for rrit = 180 
GeV, as = 0.125 and four different energies. The dotted line shows the free particle case 
Re = p/rrit (from ||12?|| ) . 



2 1 - iaJSn 
o 1 — oas/oTT 

The component of the polarization normal to the production plane is thus approximately 
independent of E and essentially measures the strong coupling constant. In fact one can 
argue that this is a unique way to get a handle on the scattering of heavy quarks through 
the QCD potential. 



Predictions for real and imaginary parts of the function ip are displayed in Fig. |3.27| for 
four different energies. 



The momentum integrated functions ^{E) are shown in Fig. |3.28| . From this figure, in 



combination with Fig. |3.26| , it is clear that the contribution of P-wave amplitudes to the 



quark polarization will amount to 10% at most and by construction vanishes upon angular 
integration. 
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Figure 3.28: Real and imaginary part of ^{E) for three different values of (from ||127|| ). 
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3.2.6 Rescattering 



For a particle with a very small decay rate production and decay amplitudes can be clearly 
separated. This is fairly evident from the space-time picture of such a sequence. Prior to 
its decay the particle travels away from the production point and any coherence is lost 
between the two reactions. The situation is different for the case under discussion, an 
unstable top quark which decays within the range of interaction between t and t. In such 
a situation the decay products from t are still affected by the force originating from t and 
vice versa (Fig. |3.29|) . 




Figure 3.29: Lowest order rescattering diagrams. 



In ref. [|128| , |129|| it has been demonstrated that the total cross section remains unaffected 



by rescattering in order a^. This result had been anticipated in ||112|| on the basis of 
earlier work which considered the decay rate of a muon bound in the strong field of a 
nucleus [ |13CI|| . In contrast momentum and angular distributions ||122| , |129| , |131|| as well as 
the top quark polarization ||127|| are affected by rescattering. For example the momentum 
distribution has to be corrected by a factor (1 + ipiip, E)) with 



MP,E)=2lmJ^^Vi\k-p\) 



G(k, E) arctan ^ 
'G(p,i^) |k-p| 



(3.91) 
(3.92) 



The distribution is shifted towards smaller momenta by about 5% (Fig. |3.3CI| ), an ef- 
fect that could become relevant in precision experiments. The influence on the forward- 
backward asymmetry and the polarization is even more pronounced ||127|| , as far as the 
S'-P-wave interference terms are concerned which are thus intrinsically of order (3. A 
detailed discussion of these effects is beyond the scope of these lectures and can be found 
in [|l27i. 



3.2.7 Relativistic corrections 

In 0{a1) one anticipates effects from relativistic corrections, from the reduction of the 
phase space through the binding energy and from the Coulomb wave function of the h 
quark. Individually these effects are large. For the sake of the argument, let us adopt 
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Figure 3.30: Modification of tlie momentum distribution tlirougli rescattering. Daslied 
line: no rescattering corrections included; Solid line: rescattering contribution with full 
potential included; dotted line: rescattering contribution with pure Coulomb potential 
and as = 0.187 included (from [|T2^ 
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a pure Coulomb potential and a binding energy of -2.5 GeV. From the virial theorem 
one derives a potential energy of -5 GeV. The phase space of the quark decaying first is 
therefore reduced by this same amount. Assuming rrit = 180 GeV one would arrive at a 
reduction of Tt by about 10%. A full calculation of all 0{a1) effects is not available at 
present and one has to resort to models and analogies ||112| , |119| , |114|| . For example, it 



has been shown ||130| , |132|| that the decay rate of a muon bound in the field of a nucleus 
is given by 



free 



5{Za) 



5{Za) 



{Zaf 



(3.93) 



where the first correction factor originates from the phase space suppression, the second 
from the Coulomb enhancement, and the third from time dilatation. Thus there is no 
first order correction to the total rate from rescattering in the nucleus potential, similar 
to the ti case discussed above. The second order contributions evidently compensate 
to a large extent. In a model calculation where these features are implemented [|114| 
through a momentum dependent width, it is found that the total cross section as well as 
the momentum distribution are hardly affected. These considerations have recently been 
confirmed in a more formal approach ||133 |. 
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